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Two  atmospheric  models  were  used  to  predict  nocturnal  cooling  in  a 
vegetated  environment.  The  results  from  both  models  were  compared  with 
the  observed  temperatures  to  determine  which  model  predicted  the  cooling 
curve  more  accurately.  Blackadar's  Boundary  Layer  Model  (ELM)  was  tested 
against  the  original  Satellite  Frost  Forecast  System  prediction  model, 
called  the  P-model.  Upper  air  soundings  obtained  from  the  National  Weather 
Service  computer  network  and  local  characteristics  of  the  nocturnal 
environment  were  processed  into  Boundary  Layer  Model  input  files.  Ground 
station  measurements  were  used  as  input  to  the  P-model,  and  remotely-sensed 
temperatures  from  Geostationary  Operational  Environmental  Satellite  (GOES 
West)  were  used  as  input  to  a  modified  P-model.  The  models  were  run  on 
a  real-time  basis  on  the  eve  of  23  advective  and  radiative  nights.  Model 
output  was  analyzed  for  two  Florida  locations,  Gainesville  and  Ruskin. 
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The  BLM  predioted  temperatures  with  more  precision  than  both  the 
?-model  and  the  modified  P-model  for  both  sites.  The  95%  confidence 
intervals  from  t-tests  run  to  determine  significant  difference  between 
predicted  minus  observed  and  zero  averaged  +  1.7  C  for  the  BLM  at  Gainesville, 
+  2.9  C  for  the  P-model  and  modified  P-model  at  Gainesville,  +  1.4  C  for 
the  BLM  at  Ruskin,  and  +  4.6  C  for  the  modified  P-model  at  Ruskin.  A 
BLM-predicted  temperature  bias  of  +  3.35  C  at  the  Gainesville  site  was 
attributed  to  the  interpolation  procedure  that  produced  a  sounding  for 
Gainesville  from  the  Waycross  and  Ruskin  soundings. 

Temperature  predictions  have  a  direct  application  to  horticultural 
users  in  frost  protection  management  decisions.  Prediction  of  temperature 
versus  time  provides  information  on  when  critical  temperatures  are  being 
reached,  the  duration  of  those  temperatures,  and  the  minimum  temperature. 
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INTRODUCTION 


Temperature  predictions  aid  in  frost  protection  decisions,  providing 
the  horticultural  user  warnings  of  critical  temperatures.  Frost  night 
temperature  forecasts  are  usually  limited  to  minimum  temperature,  but 
prediction  of  cooling  rates  in  small  time  increments  would  be  much  more 
beneficial.  A  predicted  cooling  curve  could  provide  the  user  with  not 
only  minimum  temperature,  but  also  time  of  occurrence  of  mimiraum  temperature, 
time  of  occurrence  of  critical  temperature,  and  duration  of  critical 
temperature.  The  high  cost  of  heating  oil,  wind  machine  fuel,  and  limited 
water  in  Florida  are  incentives  to  conserve  these  resources.  Readily 
available  temperature  predictions  may  enable  a  horticulturist  to  save 
hours'  worth  of  fuel  or  water  during  the  winter. 

Frost  protection  originated  thousands  of  years  ago  (Martsolf,  1979b), 
but  temperature  prediction  efforts  are  documented  in  more  recent  history 
(Georg,  1971).  The  earliest  prediction  formulae  date  back  to  the  late 
I800s.  The  incentive  to  develop  methods  of  local  temperature  predictions 
began  with  agricultural  needs,  after  severe  freeze  damage  to  agricultural 
crops  occurred  in  the  western  states.  Early  predictions  employed  simple 
empirical  formuli,  but  much  more  elaborate  theoretical  and  empirical 
methods  have  been  used  ever  since. 

The  Federal-State  Frost  Warning  Service  has  been  providing  frost 
night  temperature  predictions  on  an  operational  basis  for  over  50  years. 
The  Satellite  Frost  Forecast  System  (SFFS)  was  developed  as  an  aid  to  that 
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service  (Figure  1).  SFFS  provided  satellite-sensed  temperature  maps  of 
Florida  and  computer-produced  temperature  predictions  for  several  locations 
in  the  state  to  National  Weather  Service  forecasters  in  Ruskin,  Florida, 
where  the  frost  forecasts  were  made.  The  computer-produced  predictions 
were  calculated  from  a  model  that  simulated  the  surface  energy  balance 
under  nocturnal  conditions.  This  model,  called  the  P-model,  used  automated 
weather  station  data  as  input.  The  opportunity  arose  to  obtain  a  second 
model  capable  of  predicting  temperatures.  This  model,  Blackadar's  Boundary 
Layer  Model,  was  a  more  rigorous  treatment  of  the  atmospheric  boundary 
layer  than  the  original  SFFS  model  (Blackadar,  1979).  The  model  has  been 
used  by  Blackadar  mostly  for  educational  and  demonstration  purposes.  This 
is  an  application  of  the  model  under  a  specific  condition — cold  nights  in 
Florida. 

The  purpose  of  this  study  is  to  apply  the  Boundary  Layer  Model  to 
the  Florida  nocturnal  environment,  compare  its  predictions  and  the  P-model 
predictions  with  the  observed  cooling  curve,  and  determine  which  model 
provides  a  more  realistic  prediction  of  the  observed  cooling  curve. 
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LITERATURE  REVIEW 


This  project  was  developed  as  part  of  the  Satellite  Frost  Forecast 
System  (SFFS).  The  forecasting  aspect  of  SFFS  dealt  with  the  use  of  a 
temperature  prediction  model.  Therefore,  the  review  covers  two  areas: 
The  Satellite  Frost  Forecast  System,  to  provide  a  historical  background 
of  the  temperature  prediction  effort  in  Florida;  and  atmospheric  modeling, 
to  show  examples  of  other  applied  research  in  the  area  of  prediction  of 
atmospheric  variables. 

The  Satellite  Frost  Forecast  System 

The  Satellite  Frost  Forecast  System  was  developed  through  contracts 
with  the  National  Aeronautics  and  Space  Administration  (NASA)  ,  with  the 
original  goal  of  providing  GOES  (Geostationary  Operational  Environmental 
Satellite)  thermal  infrared  images  and  results  of  a  temperature  prediction 
model  to  National  Weather  Service  (NWS)  forecasters.  These  Images  and 
predictions  were  to  aid  the  forecasters  in  making  predictions  during  frost 
nights.  Several  reports  from  the  Climatology  Laboratory  at  the  University 
of  Florida  to  NASA  document  the  details  of  the  system  development  from 
1977  to  1983  (Martsolf,  1983a). 

Bartholic  and  Sutherland  (1978)  described  a  data  collection  station 
network.  Twelve  stations  located  around  the  state  of  Florida  collected 
hourly  measurements  of  air  temperature,  wind  speed,  soil  temperature,  dew 
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point  temperature,  and  net  radiation.  These  measurements  were  made  manually 
by  volunteers.  The  measurements  were  used  as  input  to  a  temperature 
prediction  model,  called  the  Physical  model  (or  P-model)  (Bartholic  and 
Sutherland,  1973;  Sutherland,  I98O).  The  P-model  was  developed  for  SFFS 
and  based  on  previous  efforts  by  Georg  (1971).  The  P-model  simulated  a 
soil  surface  energy  balance.  It  provided  1.5-m  air  temperature  predictions 
throughout  a  frost  night  for  several  locations  across  Florida.  Sutherland 
(1980)  reported  that  the  model  predicted  within  3  C  of  the  observed 
temperature  in  98^  of  141  runs  and  within  2  C  in  81.6J  of  the  runs  during 
the  winter  of  1977-1978.  The  use  of  synoptic  forecasts  of  radiation  and 
wind  speed  as  model  input  Increased  its  accuracy.  Further  examples  of 
P-model  results  were  given  by  Martsolf  and  Gerber  (I98I).  The  P-model 
predictions  were  considered  satisfactory,  with  mean  errors  of  -0.4  C  and 
-0.2  C  for  two  given  examples,  but  the  standard  deviations  ware  high, 
indicating  large  variability  in  the  hour-to-hour  predictions: 

The  Geostationary  Operational  Environmental  Satellite  (GOES)  East 
images  were  obtained  on  a  near  real-time  basis  to  aid  forecasters  in  frost 
temperature  prediction  (Barnett  et  al.,  I98O).  The  maps  were  originally 
obtained  from  the  National  Oceanic  and  Atmospheric  Administration  (NOAA) 
in  Miami,  Florida,  through  phone  lines  (Sutherland  et  al.,  1979;  Gaby, 
1980).  A  Hewlett-Packai-d  (HP)  minicomputer  located  at  the  NWS  forecasting 
office  in  Ruskin,  Florida,  called  Miami  hourly  and  downloaded  data  covering 
the  state  of  Florida.  The  satellite  data  were  transmitted  through  phone 
lines  by  analog  signal,  then  digitized  Into  map  form  by  the  HP  computer. 
The  data  were  sampled  In  order  to  present  an  undlstorted  Image  of  Florida, 
but  this  process  reduced  the  data  and  introduced  some  errors.  However, 
the  errors  were  insignificant  for  qualitative  evaluation  of  the  temperatures 
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over  this  large  an  area.  The  idea  of  providing  these  maps  to  agricultural 
users  was  proposed  at  this  early  date  (Bartholic  and  Sutherland,  1978). 
One  method  proposed  was  to  display  the  maps  on  a  grower's  home  television 
screen  through  use  of  electronics  similar  to  what  was  used  in  home  television 
games.    Personal  computers  were  not  yet  readily  available. 

P-model  predictions  from  ten  sites  were  used  as  input  to  a  second 
model,  called  the  S-model  (Bartholic  and  Sutherland,  1978).  The  S-model 
used  an  observed  satellite  map  as  a  base,  and  extrapolated  the  temperature 
changes  with  time  from  the  P-model  across  the  base  map  to  produce  a 
forecasted  map.  This  map  was  identical  in  appearance  to  the  observed  map, 
but  showed  predicted  instead  of  observed  temperatures  in  color-coded  ranges. 

The  acquisition  of  digital  satellite  data  from  NOAA  files  at  the 
National  Meteorological  Center  (NMC)  in  Camp  Springs,  Maryland,  began  in 
1979,  replacing  the  digital-to-analog,  analog-to-digital  procedure  (Martsolf , 
1979a).  The  reduction  of  conversions  improved  the  accuracy  of  the  incoming 
data,  although  the  maps  appeared  to  have  a  decrease  in  resolution  due  to 
the  removal  of  smoothing  between  conversions. 

Martsolf  (1980)  described  improvements  in  SFFS  from  its  early  years. 
The  collection  stations  were  automated  during  the  Fall  of  1979.  The  Ruskin 
computer  would  interrogate  a  microprocessor  on  each  automated  weather 
station  to  acquire  wind  speed  and  temperature  measurements.  The  measurements 
were  then  automatically  read  by  the  P-model. 

The  system  began  to  move  into  a  second  phase  about  this  time,  with 
the  beginning  of  dissemination  to  agricultural  users  (Martsolf,  1980). 
An  experimental  microcomputer  network  was  developed  between  two  county 
extension  offices  and  the  SFFS  computer  in  Gainesville,  allowing  the 
transfer  of  maps  to  the  agricultural  sector.     APPLE  II  microcomputers 
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called  Into  the  SFFS  HP  computer  through  modems  and  phone  lines  to  download 
the  satellite  maps.  Methods  for  direct  interface  between  agricultural 
users  and  SFFS  were  proposed. 

Satellite  data  acquisition  success  and  automated  weather  station 
access  success  were  documented  by  Martsolf  and  Gerber  (I98I).  The  success 
rate  of  acquiring  maps  from  NMC  varied  from  7%  to  100?  which  was  considered 
unsatisfactory.  Most  acquisition  failures  were  due  to  unavailability  of 
maps  in  the  queue  at  NMC.  A  direct  digital  downlink  to  GOES  at  Gainesville 
was  proposed.  The  direct  dissemination  system  to  agricultural  users 
expanded  to  six  county  extension  offices.  All  acquired  the  satellite  data 
with  APPLE  II  computers  through  phone  lines  (Martsolf,  I98I). 

The  capability  of  SFFS  map  acquisition  increased  greatly  after  the 
installation  of  a  direct  digital  downlink  with  GOES  (Martsolf,  1982). 
The  direct  link  upgraded  the  access  rate  of  maps  from  every  hour  to  every 
half  hour  and  doubled  the  number  of  pixels  per  map.  The  addition  of  the 
Institute  of  Food  and  Agricultural  Sciences  (IFAS)  VAX  750  computer  as 
another  node  in  the  system  was  also  described  by  Mairtsolf .  The  information 
was  transmitted  to  the  VAX  on  a  regular  basis  through  a  4800-baud  dedicated 
line.  The  VAX  allowed  anyone  associated  with  IFAS  that  had  access  to  a 
microcomputer  or  terminal  to  obtain  SFFS  weather  data.  The  number  of 
agricultural  extension  offices  in  the  network  grew  to  10.  Martsolf  mentions 
the  possible  use  of  a  Control  Data  Corporation  CYBER  730  to  aid  in 
acquisition  and  dissemination. 

Jackson  and  Ferguson  (1983)  developed  an  experimental  microcomputer 
network  to  disseminate  SFFS  information  and  other  products  to  growers. 
The  system,  called  LOIS  (Lake-Orange  Information  System),  used  the  county 
extension  office  APPLE  11+  microcomputers  as  nodes  to  growers  who  had 
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their  own  microcomputers.  Initially,  eight  growers  would  call  the  extension 
offices  through  modems  and  phone  lines  to  download  relevant  information. 
Users  could  call  in  as  often  as  they  wished  except  during  freeze  nights 
when  demand  increased  considerably. 

The  Satellite  Frost  Forecast  System  began  acquiring  text  weather 
information  from  NWS  through  a  dedicated  phone  line  between  Ruskin  and 
the  SFFS  computer  In  Gainesville  (Martsolf,  1983a).  A  listen-only  link 
was  made  between  the  NWS  HP  computer  in  Ruskin  and  the  Automated  Field 
Operational  System  (AFOS),  which  is  the  NWS  computer  network  through  which 
weather  information  is  passed.  This  link  was  initially  established  to 
obtain  dew  point  temperatures  from  the  State  Weather  Roundups  to  be  used 
as  an  additional  input  to  the  P-model  (Martsolf,  1983b).  It  was  found 
that  the  dew  points  were  still  sufficiently  high  in  the  early  evening  of 
a  frost  night  to  cause  the  P-model  to  underpredict.  The  use  of  forecasted 
dewpoint  minima  was  tried  with  more  success.  The  NWS  fruit  frost  forecasters 
requested  that  the  system  permit  the  delivery  of  their  minimum  temperature 
forecast  from  AFOS  to  SFFS  users  through  this  link.  The  establishment  of 
the  AFOS  link  also  allowed  other  text  information  to  be  obtained  via  the 
same  link.  The  text  data  included  zone,  agricultural  weather  and  freeze 
forecasts,  weather  roundups,  and  radar  summaries. 

An  experiment  to  use  satellite  data  as  input  to  the  P-model  was  made 
by  Heinemann  and  Martsolf  (1984).  The  average  of  nine  pixels  from  each 
automated  weather  station  location  was  calculated  and  used  as  surface 
temperature.  The  success  of  the  predictions  was  limited;  the  procedure 
worked  for  areas  that  were  clear  and  dry,  but  the  presence  of  clouds  or 
moisture  interfered  with  the  correct  measurements  of  surface  temperatures, 
reducing  the  prediction  accuracy. 
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Development  of  rainfall  products  began  on  the  system  (Martsolf  et 
al.,  1984;  Heinemann  et  al.,  1984).  Manually  digitized  radar  (MDR)  data 
from  seven  radar  stations  covering  the  state  of  Florida  were  acquired 
through  AFOS.  The  cold  cloud  top  temperatures  shown  by  the  infrared 
satellite  measuranents  usually  indicated  an  area  of  rain;  the  low  resolution 
MDR  (32  km)  verified  the  presence  and  intensity  of  rain.  Overlay  of  the 
satellite  data  onto  the  MDR  increased  the  rainfall  resolution  to  6  x  8  km 
areas . 

Martsolf  (1983b)  gave  a  description  of  the  possible  successor  to 
SFFS,  Florida  Agricultural  Services  and  Technology,  Inc.  (FAST),  which 
was  designed  to  be  a  not-for-profit  organization  to  disseminate  weather 
and  other  products  to  users  not  associated  with  IFAS.  The  development  of 
products  such  as  the  rainfall  maps  and  acquisition  of  the  text  information 
made  this  effort  feasible.  They  began  using  a  CYBER  730  to  act  as  a  node 
to  the  public. 

Florida  Agricultural  Services  and  Technology  was  in  operation  by  the 
Spring  of  1984  (Martsolf,  et  al.  1985).  In  the  near  future,  the  FAST 
CYBER  is  expected  to  become  the  hub  of  the  weather  information  acquisition 
and  dissemination  system,  which  will  be  called  FAWN  (Florida  Agricultural 
Weather  Network). 

Modeling 

Modeling  is  a  useful  method  of  gaining  a  better  understanding  of  real 
world  processes.  This  review  covers  models  that  pertain  to  the  atmospheric 
boundary  layer  processes,  temperature  prediction,  and  large  scale  atmospheric 
predictions. 
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Atmospheric  prediction  models  can  be  divided  into  several  categories. 

The  atmospheric  processes  affecting  global  circulation  actually  begin  at 

the  microscale  (less  than  1  km)  and  extend  through  the  synoptic  scale 

(greater  than  100  km).    Kolsky  (1972)  estimated  that  a  model  that  encompasses 

all  aspects  of  atmospheric  processes  would  need  a  scale  ratio  on  the  order 
10 

of  10  ,  which  even  with  today's  supercomputers  would  be  unreasonable. 
Therefore,  the  physical  models  are  divided  into  synoptic,  mesoscale,  and 
microscale  regimes. 

The  first  major  large  scale  numerical  weather  forecasts  were  attempted 
by  Richardson  (1922).  Primitive  hydrodynaraic  equations  formed  the  base 
of  his  model.  Richardson  proposed  applying  the  model  to  an  operational 
weather  network,  but  there  were  two  major  problems  with  Richardson's 
approach.  First,  poor  data  for  the  initial  conditions  and  inclusion  of 
gravity  and  sound  waves,  referred  to  as  'noise'  by  Charney  (19^9),  greatly 
amplified  errors  through  the  forecast  period.  Second,  in  that  pre-coraputer 
era,  the  numerical  process  was  quite  cumbersome,  so  no  further  serious 
efforts  were  made  for  another  twenty  years  (Petterssen,  1956). 

The  invention  of  computers  made  simulations  of  physical  processes 
feasible  and  practical.  Complex  models  have  been  developed  that  can  be 
run  in  reasonable  time  periods  with  the  increased  speed  and  size  of 
computers. 

Charney  (19^9)  was  responsible  for  filtering  out  meteorological  noise 
by  reducing  the  complexity  of  the  equations.  He  determined  that  the 
high-speed  gravity  and  sound  waves  had  no  significant  effect  on  the  larger 
and  slower  moving  atmospheric  motions.  The  first  operational  weather 
prediction  model  was  the  result  of  this  work,  and  was  known  as  the  barotropic 
model.     The  barotropic  model  was  a  dynamic  model  based  on  continuity 
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equations,  but  all  pressure  surfaces  were  parallel;  density,  temperature, 
and  pressure  coinoided.  Because  of  this,  the  barotropic  model  neglected 
thermal  advection,  so  new  weather  systems  could  not  be  developed  by  the 
model . 

The  limitations  in  the  barotropic  model  led  to  the  development  of 
the  baroclinic  model  (Haltiner,  1971).  The  baroclinic  model  contained 
more  than  one  level  to  determine  thicknesses  between  pressure  surfaces, 
and  the  wind  flow  could  cross  isobars.  This  enabled  the  model  to  account 
for  temperature  advection  and  the  development  of  new  weather  systems. 

Combinations  of  these  models  were  used  to  produce  operational  weather 
forecasts.  Two  operational  weather  prediction  models  are  used  by  the 
National  Meteorological  Center.  A  six-layer  baroclinic  model  based  on 
primitive  hydrostatic,  hydrodynamic ,  and  thermodynamic  equations  was 
developed  by  Shuman  and  Hovermale  (1968).  The  model  gridding  was  380  km 
over  the  Northern  Hemisphere.  This  model  significantly  improved  the 
large-scale  weather  predictions  over  previous  efforts.  A  finer  grid  model 
(190  km)  known  as  the  Limited  Fine  Mesh  (LFM)  is  used  over  the  United 
States  and  adjacent  coastal  waters  (Holton,  1979). 

A  more  complex  forecasting  model  became  operational  in  Europe,  called 
the  SCMWF  (European  Centre  for  Medium  Range  Weather  Forecasts)  model 
(Tiedtke,  1983).  The  ECMWF  model  consists  of  15  levels.  The  grid  spacing 
is  presently  170  km,  which  is  being  improved  to  100  km  with  the  acquisition 
of  a  Cray  XMP-22  (Kerr,  1985b).  The  synoptic  model  also  includes  mesoscale 
and  microscale  processes  such  as  radiation,  condensation,  cumulus  convection, 
and  turbulence. 

The  latest  model  developed  by  the  NMC  is  the  Regional  Analysis  Forecast 
System  (RAFS)  (Kerr,  1985a),  which  has  16  vertical  layers,  and  uses  nested 
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grid  spacing.  The  Western  Hemisphere  above  the  equator  has  a  grid  spacing 
of  320  km.  A  large  square  of  160  km  is  embedded  within  the  320-km  spacing, 
and  a  second  embedded  square  of  80-km  grid  spacing  covers  North  America. 
The  RAFS  model  improved  prediction  of  some  important  processes  over  the 
LFM.  Two  examples  are  prediction  of  winter  rainfall,  both  in  coverage 
and  total  amounts,  and  the  January  1985  freeze  in  Florida. 

The  large  scale  prediction  models  generally  predict  events  occurring 
over  thousands  of  kilometers  spatially  and  several  hours  to  several  days 
temporally.  Therefore,  the  output  from  these  models  feature  synoptic 
processes  such  as  movement  and  development  of  pressure  systems.  Output 
that  predicts  some  of  the  smaller  scale  phenomena  tends  to  be  suppressed 
because  of  the  volume  of  data  produced. 

Mesoscale  (1  km  to  100  km)  and  microscale  models  emphasize  the  smaller 
scale  processes.  Temporally,  the  models  usually  do  not  predict  beyond  2k 
hours.  The  numerical  procedure  time  increments  are  in  minutes  or  seconds. 
Proper  integration  of  the  smaller  scale  models  into  the  synoptic  scale 
models  can  increase  the  accuracy  of  the  large-scale  predictions,  but  the 
limiting  factor  is  computer  capability.  The  following  paragraphs  provide 
some  examples  of  mesoscale  and  microscale  models. 

One  of  the  earlier  mesoscale  models  was  developed  by  Pielke  (197^). 
Pielke  modeled  the  sea  breeze  to  predict  shower  activity  along  the 
convergence  regions  of  Florida  coastal  winds.  It  was  based  on  a  moisture 
continuity  equation  and  the  horizontal  and  vertical  equations  of  motion 
similar  to  those  used  in  the  synoptic  scale  models,  but  Pielke 's  model 
required  greater  resolution  than  that  used  in  the  large  scale  models. 
Some  of  the  boundary  layer  parameterizations  used  were  developed  by 
Blackadar  and  Tennekes  (1958). 
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The  SFFS  P-aodel  was  based  on  a  surface  energy  balance  model  developed 
by  Georg  (1971).  The  Georg  model  simulated  the  balance  between  soil  heat 
flux,  surface  radiation  loss,  background  radiation,  and  convective  heat 
exchange  in  the  air  layer  just  above  the  surface.  The  model  used  inputs 
from  four  soil  temperatures,  a  surface  temperature,  a  1.5  m  air  temperature, 
wind  speeds  at  9  and  18  m,  net  radiation,  and  dew  point  temperature.  The 
model  produced  good  minimum  temperature  predictions  when  compared  to 
observed  values;  twenty-six  of  the  45  predicted  minimums  were  within  1  C 
of  the  observed. 

Deardorff  (1977)  suggested  a  method  of  estimating  the  soil  moisture 
fraction  (measured  soil  moisture/saturation  value)  based  on  a  rate  equation. 
Most  models  involving  soil  properties  usually  hold  the  soil  moisture 
constant.  This  method  gave  the  diurnal  variation  of  soil  moisture  fraction 
over  several  days.  It  involved  bulk  soil  moisture  content,  soil-surface 
evaporation  rate,  and  precipitation  rate. 

A  one-dimensional  numerical  model  simulating  surface  temperature  and 
heat  flux  was  developed  by  Carlson  and  Boland  (1978).  This  model  used 
K-type  parameterization,  included  radiative  flux  divergence  (0.15  C  h  ^ 
at  night),  and  neglected  advection.  The  model  describes  a  1.5-meter  deep 
substrate  beneath  the  earth's  surface,  a  50-meter  surface  layer  which 
includes  a  transition  layer  over  a  crop  canopy  or  urban  surface,  and  a 
mixed  layer  above.  They  concluded  that  the  thermal  inertia  and  the  moisture 
availability  were  the  two  most  important  parameters  in  this  model,  but 
that  these  are  difficult  to  measure  for  heterogeneous  urban  and  rural 
terrain. 

A  forest  microclimate  was  modeled  by  Waggoner  et  al.  (1969).  The 
forest  canopy  energy  exchange  was  simulated,  based  on  radiant,  sensible, 
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and  latent  heat  fluxes.  The  model  was  used  to  study  effects  of  stomatal 
changes  on  evaporation  within  the  forest  canopy.  This  model  used  an 
electrical  network  analogy,  representing  the  various  fluxes  as  flow  and 
resistances. 

Some  prediction  models  were  based  on  statistical  rather  than  physical 
considerations.  Aron  (1975)  used  regression  equations  to  predict  the 
number  of  hours  that  air  temperature  would  remain  below  a  given  threshold, 
in  order  to  calculate  chilling  units.  The  regressions  were  based  on  a 
16-year  sample  of  hourly  temperature  and  humidity  measurements  recorded 
at  first  order  collection  stations  in  California.  Linear  regression  of 
predicted  versus  observed  hours  below  the  threshold  temperature  of  45  F 
gave  a  correlation  coefficient  on  the  order  of  0.98. 

Petersen  (1976)  described  a  model  that  simulates  atmospheric  turbulence,  • 
through  the  use  of  generalized  spectral  analysis.    This  model  was  based 
on  a  statistical  approach.    The  model  simultaneously  generated  all  three 
velocity  components,  and  dealt  with  velocity  variations  on  the  order  of 
seconds . 

A  model  to  predict  hourly  temperatures  through  a  purely  statistical 
method  was  developed  by  Hansen  and  Driscoll  (1977).  The  parameters  are 
based  on  an  11-year  developmental  sample.  Hansen  and  Driscoll  concluded 
that  the  model  is  most  effective  for  predictions  of  long-term  tenperature 
events,  such  as  durations  of  temperatures  above  or  below  a  certain  point, 
and  temperature  extremes. 

Parton  and  Logan  (1981)  developed  a  model  that  predicted  diurnal 
variation  in  air  and  soil  temperatures.  The  model  used  mathematical 
functions  to  determine  temperature  based  on  maxijna  and  minima.  A  truncated 
sine  wave  was  used  for  daytime  temperatures  and  an  exponential  function 
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was  used  for  nighttime  temperatures.    The  results  produced  a  maximum 
absolute  mean  error  of  2.64  C  for  a  10-cm  air  temperature. 

Another  statistical  temperature  prediction  model  was  developed  by 
Miller  (1981).  This  model,  known  as  GEM,  predicted  temperatures  and  other 
meteorological  variables  up  to  12  hours  in  advance.  The  predictions  were 
calculated  through  multiple  linear  regression  of  local  meteorological 
observations.  The  GEM  model  calculated  the  probability  that  a  temperature 
will  fall  within  a  certain  range  for  each  hour. 


MODELS  AND  METHODS 

The  two  models  used  in  this  study  were  the  P-model,  which  was  the 
original  SFFS  prediction  model,  and  Blackadar's  Boundary  Layer  Model  (BLM). 
The  models  were  used  to  predict  a  cooling  curve  of  a  nocturnal  environment. 
Although  the  two  models  are  of  different  scales  (BLM  is  mesoscale,  the 
P-model  is  microscale),  the  energy  balance  approach  used  in  the  P-model 
is  similar  to  some  of  the  low-level  microscale  processes  simulated  by  the 
BLM.  This  and  the  willingness  and  cooperation  of  the  BLM's  author  in 
allowing  the  BLM  to  be  used  for  this  study  were  the  main  reasons  for  the 
selection  of  this  particular  model.  The  BLM  was  obtained  directly  from 
Dr.  Alfred  Blackadar  by  coraputer-to-coraputer  link  over  phone  lines 
(Blackadar,  1983>  personal  communication).  This  section  describes  the 
important  physics  and  equations  of  each  model.  Each  model  is  divided  into 
discernible  modules,  each  module  representing  a  physical  process  simulated 
by  the  model. 

The  Models 


The  P-model 

The  P-model  sinulates  a  soil  surface  energy  balance  (Sutherland, 
1980)  (Figure  2).  The  model  is  based  on  previous  work  by  Georg  (1971). 
The  model  can  be  divided  into  three  modules,  each  representing  the  rate 
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Figure  2.    Physical  processes  involved  in  the  P-model. 
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of  temperature  change  with  time,  and  each  is  governed  by  a  differential 
equation: 

the   3ub-3urface  layer, 

s  at       33z^  ' 

the   earth-air  interface, 

C„   ~  =       +   S   +   R   +   B    ,  (2) 

s   d  t 

and   the   air  layer, 

3_T  H  1        dU»  H  

3t   "   "   h         aU,     dt     K(U,,,z)  ' 

where 

-1  -1 

C    =  soil  specific  heat  (KJ  Kg      C    ) , 

3 

T    =  soil  temperature  (C), 

3 

t  =  time  (s), 

z  =  depth  or  height  (m), 

-2 

H  =  air  layer  sensible  heat  flux  (W  m    ) , 

_2 

S  =  soil  heat  flux  (W  m    ) , 

_2 

R  =  surface  radiative  flux  (W  m    ) , 

_2 

B  =  background  radiative  flux  (W  m    ) , 

1  ? 

C    =  air  volumetric  specific  heat  (J  K     m  ), 

h  =  parameter  related  to  height  of  nocturnal  boundary  layer  (m) , 
U„  =  friction  velocity  (m  s~^), 

=  convective  heat  transfer  coefficient  (W  m~    C~  ) , 

1  1 

K    =  thermal  conductivity  (W  m     C~  ) . 


Sutherland  does  not  define  a ;  this  is  discussed  in  following 
paragraphs. 
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Equation  1  is  the  Fourier  heat  transfer  equation  for  conduction  of 
energy  through  a  solid.  Equation  1  is  applied  to  the  soil  in  the  P-model; 
numerical  approximations  of  the  Fourier  equation  are  used  to  determine 
the  soil  heat  flux  from  a  depth  of  constant  temperature  to  the  surface 
(Wierenga  and  de  Wit,  1970). 

Equation  2  is  the  surface  energy  balance,  which  drives  the  model. 

This  form  of  the  energy  balance  is  non-steady  state  but  would  approach 

steady  state  as  the  temperature  change  with  time  term  on  the  left-hand 

side  approaches  zero.    The  soil  heat  flux  S  is  provided  by  equation  1. 

The  air  layer  sensible  heat  flux  found  in  equation  2  is  the  same  value 

calculated  from  equation  3.  The  surface  radiation  R  is  calculated  from 

the  Stephan-Boltzmann  law  using  the  ground  surface  temperature  T  ,  and 

3  ' 

the  background  radiation  B  is  calculated  from  initial  measurements  of  net 
radiation. 

The  sensible  heat  flux  H  in  equation  3  is  calculated  from 


AT 

H   =    K(U,,z)--  , 
Az 

K(U,,z)    =   K^expCSU.z)  , 


(4) 


(5) 


and 


U,   =  ka-ln(z/z^) 


(6) 


where 


Ug  =  wind  speed  at  9m  (m  s~  ) 
z  =  height  (m), 


z^  =  roughness  length  (m), 


k  =  Von  Karman  Constant  (0.4). 


Beta  and  are  derived  by  setting  the  right-hand  side  of  equation  5  equal 
to  the  linear  form  K  =  kU,z. 
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There  was  controversy  with  regard  to  Sutherland's  formulation  of 

some  model  components.    Shaw  (1981)  commented  on  four  sections  of  the 

P-model.    First,  Shaw  questioned  the  possibility  of  an  infinitely  thin 

surface  having  a  mass,  as  characterized  by  C    in  equation  2.     The  units 

s 

of  C    were  also  inconsistent;  Shaw  contends  that  the  heat  capacity  per 
s 

unit  area  for  the  surface  should  not  have  the  same  symbol  (and  therefore 

the  same  units)  as  the  soil  heat  capacity  per  unit  volume.  Sutherland 

(1981)  counters  that  equation  2  is  in  steady-state,  therefore  the  left 

side  becomes  zero  and  the  solution  is  trivial.     Sutherland  could  have 

avoided  this  problem  in  two  ways.    The  surface  would  have  a  volumetric 

heat  capacity  if  a  finite  surface  thickness  of  perhaps  a  few  millimeters 

was  used.    The  steady  state  assumption  would  still  set  the  left  side  to 

zero,  so  the  presence  of  a  volumetric  heat  capacity  would  not  change  the 

outcome  of  the  energy  balance  calculations.    Also,  by  using  a  flux  form 

of  the  left  side  of  equation  2,  the  need  for  a  heat  capacity  term  would 

have  been  eliminated. 

Second,  a  question  arises  from  Sutherland's  wind  profile  equation 

(right-hand  side  of  equation  3),  since  the  parameter  a  is  not  defined  for 

this  equation.    Sutherland  points  out  that  this  is  a  typographical  error, 

2 

and  the  term  should  be  gU,j  instead  of  aU,. 

Third,  a  related  question  is  the  use  of  the  exponential  wind  profile 
that  had  been  derived  for  a  plant  canopy  when  no  mention  of  a  plant  canopy 
was  made.  Sutherland  justifies  its  use  because  it  works  better  than  the 
neutral  form,  but  he  would  have  been  more  convincing  if  he  had  mentioned 
that  the  P-model  was  being  applied  to  a  vegetated  environment.  Another 
problem  with  the  development  of  the  convective  heat  transfer  coefficient 
equations  is  that  Sutherland  uses  a  linear  wind  profile  form  in  equation 
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5  to  determine  the  parametei^s  used  In  the  exponential  wind  profile  (equation 
3).  The  P-model  is  modele  use  only  a  one-level  wind  input  and  gin  assumed 
roughness  height  where  the  wind  velocity  should  be  zero.  Two  measured 
wind  levels  and  the  roughness  height  would  provide  a  better  value  of  the 
oonvective  heat  transfer  coefficient. 

A  fourth  and  final  point  of  contention  is  the  use  of  the  surface 
temperatures  for  determination  of  soil  and  sensible  heat  fluxes  since  the 
near-surface  temperature  profile  would  be  highly  non-linear.  This  is  true 
within  a  few  centimeters  of  the  surface,  but  measurement  of  the  temperature 
profile  within  a  few  centimeters  of  the  surface  is  very  difficult.  More 
importantly,  the  error  imposed  from  the  assumption  that  the  surface  and 
the  air  just  above  the  surface  are  at  the  same  temperature  would  probably 
not  be  significant  enough  to  cause  a  significant  error  in  the  flux 
calculations. 

« 

The  Modified  P-model 

The  original  version  of  the  P-model  used  automated  weather  station 
(AWS)  data  as  input.  Ten  stations  were  located  around  the  state  of  Florida. 
The  SFFS  computer  would  call  the  stations  each  hour  to  acquire  the  input 
measurements.  Two  major  problems  arose  from  the  automated  weather  station 
data  acquisition  that  reduced  the  data  reliability  (Martsolf,  1983a). 
First,  high  noise  levels  on  the  phone  lines  between  the  AWS  and  the  main 
computer  interfered  with  the  received  signal  quality.  Second,  the  AWS 
hardware  sometimes  failed.  The  hardware  failure  included  the  AWS 
microprocessor,  modem,  or  measurement  sensors.  The  AWS  data  acquisition 
success  rate  was  approximately  70%  during  the  last  winter  of  full  AWS 
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operation.  The  idea  of  using  remotely-sensed  temperatures  from  satellite 
as  input  to  the  P-model  was  proposed  as  an  alternative  to  the  automated 
weather  station  measurements  (Heinemann  and  Martsolf,  1984). 

The  P-model  was  modified  to  accept  GOES  (Geostationary  Operational 
Environmental  Satellite)  surface  temperatures  as  input  instead  of  AWS 
measurements.  The  model  ran  correctly  only  if  the  AWS  inputs  indicated 
the  presence  of  a  nocturnal  inversion.  Restrictions  included  in  the 
previous  P-model  developed  an  inversion  in  the  event  that  automated  weather 
station  data  were  missing  or  were  insubstantial.  A  subroutine  within  the 
model  read  the  surface  temperature,  and  added  1.0  C  for  the  1.5  m  temperature, 
2.0  C  for  the  3  m  temperature,  and  3.0  C  for  the  9.0  m  temperature,  but 
only  if  the  measured  change  in  temperature  with  height  was  less  than  this. 
In  the  case  of  the  Modified  P-model,  the  model  read  the  satellite  surface 
temperature,  then  used  these  restrictions  to  produce  an  inversion  profile 
based  on  the  surface  temperature. 

The  Boundary  Layer  Model 

Blackadar's  Boundary  Layer  Model  (BLM)  simulates  the  boundary  layer 
processes  occurring  in  the  first  2000  meters  of  the  atmosphere  over  a 
24-hour  period  (Blackadar,  1976).  Although  the  BLM  is  a  mesoscale  model, 
it  rigorously  models  many  of  the  microscale  physical  processes  necessary 
for  prediction  of  the  low-level  temperatures.  The  BLM  is  built  upon 
theoretical  and  empirical  relationships,  and  is  based  on  equations  for 
the  mean  velocity  (u  and  v  components),  and  the  mean  potential  temperature 
at  each  of  30  layers  (assuming  vertical  velocity  is  zero): 

St  3x  9y  3z  3x  3z      m9z    '  ^'^ 
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9V       „3V  ,,9V 

- —  +  U— —  +  V——  +  w 

d  t  9  X  9  y 


(8) 


90       „90       „90       „9  0  19R  9^90 

+   u--  +   Vr--  +  W--  =  -n  +  K. 

9t         9x         9y         9z         Cp9z  9zh3z 


(9) 


where 


-1 


U,V  =  east-west  and  north-south  wind  components  (m  s~  ), 
t  =  time  (s), 

x,y,z  =  spatial  coordinates  (ra) , 
f  =  coriolis  parameter  (s"''), 

2  -1 

=  turbulent  exchange  coefficient  for  momentum  (m    s    ) , 
p  =  air  density  (kg  m  ), 

2    - 1 

=  turbulent  exchange  coefficient  for  heat  (m    s    ) , 

-1  -3 

Cp  =  air  heat  capacity  (kJ  C     m    ) , 

0  =  potential  temperature  (K), 

_2 

R    =  Net  Radiation  (W  m    ) , 
n 

_2 

p  =  atmospheric  pressure  (n  m  ). 

Equations  7  and  8  describe  the  momentum  transfer  processes.  Equation 
7  is  the  east-west  component  and  equation  8  is  the  north-south  component. 
The  first  term  on  the  left  side  of  each  of  these  two  equations  is  the 
change  in  wind  speed  with  change  in  time.  The  second,  third,  and  fourth 
terms  on  the  left  side  of  these  two  equations  are  the  mass  advection  terms. 
The  first  term  on  the  right  side  is  the  pressure  gradient.  The  laat  terra 
on  the  right  side  describes  the  vertical  turbulent  eddy  momentum  transfer. 

Equation  9  describes  the  heat  transfer  processes.  The  first  term  on 
the  left  side  of  equation  9  is  the  rate  of  temperature  change  with  time. 
The  surface  air  layer  temperature  change  with  time  is  the  desired  output 
for  this  study.    The  second,  third,  and  fourth  terms  on  the  left  side  are 
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the  advection  terms.  The  first  terra  on  the  right  side  of  equation  9 
accounts  for  the  radiative  transfer.  The  last  term  on  the  right  side 
describes  the  vertical  turbulent  eddy  heat  transfer. 

The  BLM  can  be  divided  into  modules,  each  providing  an  essential 
process  to  the  model.  A  diagram  of  the  BLM  modules  is  shown  in  Figure  3. 
The  modules  are  integrated  so  that  processes  occuring  in  one  may  have 
direct  effect  on  another.  Blackadar  (1979)  provides  a  detailed  discussion 
of  the  model. 

The  main  purpose  of  the  soil  slab  model  is  to  provide  a  soil  surface 

temperature,  T^,  for  calculation  of  radiative  and  sensible  heat  fluxes 

from  the  surface  (Figure  4).    Conduction  through  the  soil  is  calculated 

from  a  solution  of  the  Fourier  heat  transfer  equation  (equation  1),  under 

the  conditions  that  at  an  infinite  depth  the  soil  temperature  approaches 

a  constant  temperature  T^  and  that  the  soil  heat  flux  is  continuous  at 

the  surface.    The  conductive  flux  through  the  soil  is  dependant  on  soil 

conductivity,  thermal  diffusivity,  heat  capacity,  and  water  content. 

T    is  calculated  from 
g 

~g  =  --(S  +  l4--aT'*-H)-i.<(T  -  T),  (10) 
3tC  go  mgm'' 

g 

where 

_2 

S  =  solar  insolation  (W  m    ) , 

_2 

1+  =  long  wave  atmospheric  back  radiation  (W  m  ), 

-2  -4 

a  =  Stephan-Boltzmann  constant  (W  m     K    ) , 

=  heat  flux  removed  from  surface  by  turbulence  (W  m~  ) , 

'<    =  1.18  , 
m  ' 

(li  =  angular  velocity  of  earth's  rotation  (s"^), 

-1  -3 

=  soil  heat  capacity  (kJ  C      m    ) . 
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Figure  3.     The  Boundary  Layer  Model  modules  and  important  aspects 
of  each.     T    is  the  surface  air  layer  temperature. 
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Figure  4.     The  BLM  soil  slab  model.     T    and  T    are  from  equation  10. 
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The  vegetative  model  accounts  for  water  transport  from  the  soil  to 
the  air  by  way  of  the  plant,  and  is  affected  by  root,  stomatal,  and  leaf 
boundary  layer  aerodynamic  resistances  (Figure  5).  Vegetative  radiation 
and  sensible  heat  fluxes  are  also  calculated.  Removal  of  water  from  the 
soil  affects  the  soil  heat  flux  by  lowering  the  total  soil  heat  capacity. 

Daily  radiation  loss  and  large-scale  subsidence  are  calculated.  The 
radiation  loss  would  occur  through  the  fall  and  early  winter,  until  the 
radiation  balance  of  the  northern  hemisphere  reverses  and  net  heating 
begins.  The  subsidence  term  is  added  as  a  cooling  mechanism.  After  the 
passage  of  a  cold  front,  air  is  sinking  as  the  high  pressure  builds  in 
behind  it.  The  sinking  air  would  experience  adiabatic  warming,  but  the 
net  effect  of  the  subsiding  air  on  the  lower  levels  of  the  boundary  layer 
would  be  cooling.  These  are  simply  additive  terms,  so  a  constant  value 
is  provided  for  each  process. 

Turbulent  exchange  accounts  for  the  heat,  moisture,  and  momentum 

transfer  between  the  upper  air  layers  (Figure  6).    Turbulent  transfer  is 

described  by  K-type  parameterization.     The  turbulent  heat  transfer 

coefficient  (K^)  and  the  coefficient  of  momentum  transfer  (K  )  are  assumed 
h  ni  - 

equal.    The  K  values  are  calculated  from  the  equation 


Rc  - 


Ri 


if   Ri  >    Rc , 


(11) 


K 


=  1.1 


m 


Rc 


=    0.0      if   Ri    <  Rc, 


where 


Richardson  number, 


He 


0.25  (critical  Richardson  number), 


1 


parameter  that  characterizes  the  size  of  turbulent  eddy  (m) , 


s 


wind  shear  (m  s"  /m) . 
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Figure  5.     The  BLM  vegetation  model. 
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UPPER  AIR  LAYERS 
(29  TOTAL) 


— 


SURFACE  AIR  LAYER 


Figure  6.  Physical  representation  of  the  upper  air  layers,  showing 
turbulent  heat,  momentum,  and  moisture  exchange.  Km  and 
Kh  are  from  equation  11;  terms  using  T,  U,  and  Q  are  the 
heat,  wind,  and  moisture  gradients,  respectively. 
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The  Richardson  number  is  the  ratio  of  bouyant  to  mechanical  forces, 
and  is  used  as  an  indicator  of  atmospheric  stability.  It  includes  wind 
speed,  mixing  length,  and  temperature  gradient  (Rosenberg,  1974;  Monteith, 
1980).  Generally,  the  Richardson  number  distinguishes  between  stability 
and  instability,  and  dominance  of  mechanical  or  bouyant  turbulence  within 
the  unstable  regime.  The  sign  of  the  Richardson  number  is  determined  by 
the  temperature  gradient;  positive  indicates  inversion  conditions  and 
negative  indicates  lapse  conditions.  Large  negative  values  characterize 
strong  instability  due  to  predominantly  bouyant  forces,  small  positive 
values  (less  than  0.25  in  this  model)  characterize  dominance  by  mechanical 
turbulence,  large  positive  values  characterize  strong  stability  and 
suppressed  turbulence.  The  atmosphere  in  the  boundary  layer  is  usually 
under  a  stable  turbulent  regime  during  clear  calm  nocturnal  conditions, 
with  Ri  just  less  than  0.25. 

The  heat  transfer  between  the  surface  air  layer  and  the  first  upper 
air  layer  above  is  turbulent  (Figure  7).  However,  heat  transfer  in  the 
surface  air  layer  is  complicated  by  other  transfer  processes  such  as 
sensible  heat  resulting  from  soil  heat  flux,  radiation  loss  from  the 
surface,  moisture  from  soil  and  vegetation,  and  radiation  and  sensible 
heat  from  vegetation. 

The  model  uses  an  assumed  height  of  1.0  meter  above  which  temperature 
change  is  dominated  by  turbulence  and  below  which  is  dominated  by  radiative 
flux  divergence  and  turbulent  flux  convergence.  The  surface  air  layer 
sensible  heat  flux  is  calculated  from  determinations  of  T,  and  u,,  under 
normal  stable  turbulent  conditions  occurring  at  night: 
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Figure  7.    Physical  representation  of  processes  occurring  within 
the  surface  air  layer. 
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kW 


u 


ln(z,/z 
1  o 


(13) 


m 


=  5z,/L 


(14) 


where 

=  temperature  in  first  upper  air  layer  (C), 
k  =  von  Karman  constant  (0.4), 

=  wind  velocity  in  first  upper  air  layer  (m  s"^), 

=  roughness  length  (m) , 
z^  =  height  of  surface  air  layer  (m). 

L  is  the  Monin  length,  which  is  a  function  of  the  heat  and  momentum  fluxes. 
As  with  the  Richardson  number,  it  is  a  parameter  that  indicates  stability 
by  measuring  the  ratio  of  energy  produced  by  buoyant  forces  and  energy 
dissipated  through  mechanical  forces  (Monteith,  1980). 

The  BLM  was  adapted  to  the  nocturnal  environment  (Heinemann  et  al., 
1985).  The  one-dimensional  version  used  in  this  study  neglected  the  mass 
advection  terms  found  on  the  left-hand-side  of  equations  7  and  3,  and  the 
heat  advection  term  found  in  equation  9-  The  model  contains  a  free 
convection  section  that  was  not  considered  in  this  project,  since  under 
frost  conditions  free  convection  is  very  unlikely  to  occur.  Free  convection 
occurs  when  strong  surface  heating  creates  positive  buoyant  forces  in  the 
surface  air  layer,  forcing  large  parcels  of  air  to  ascend  into  the  upper 
layers.  Under  radiative  frost  conditions,  stratification  of  the  lower 
boundary  layer  and  the  development  of  the  nocturnal  inversion  cause  the 
buoyancy  forces  to  become  negative,  hence  an  increase  in  the  lower  level 
stability  (Monteith,  1980),  but  turbulence  from  mechanical  forces  is  still 
prevalent. 
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Model  Inputs 

Each  model  requires  a  set  of  initial  Gonditions  provided  by  an  input 
file.  The  P-model  and  modified  P-model  use  three  sets  of  conditions 
observed  hourly  as  their  prediction  base.  The  BLM  uses  one  set  of  conditions 
as  input. 

P-model  Inputs 

The  inputs  required  to  drive  the  P-model  are  shown  in  Figure  8.  Soil 
heat  flux  was  determined  from  the  differential  in  the  initial  soil 
temperatures.  Outgoing  radiation  was  calculated  from  the  measured  surface 
temperature  through  the  Stephan-Boltznann  equation,  and  background  radiation 
was  measured  by  net  radiometer.  The  air  temperatures  provide  a  profile 
of  the  lower  section  of  the  inversion  layer.  The  convective  heat  exchange 
was  calculated  from  the  initial  wind  speed  and  air  temperature  profile. 

Modified  P-model  Inputs 

The  satellite  data  used  in  this  study  were  obtained  from  GOES  West 
meteorological  satellite  for  the  winter  of  1984-85  (after  the  failure  of 
GOES  East  in  July  1984).  The  Geostationary  Operational  Environmental 
Satellite  is  positioned  in  geosynchronous  orbit  38,000  km  above  the  equator 
over  the  western  hemisphere  (Anonymous,  1980).  The  satellite  is  centered 
over  the  equator  and  longitude  100  W,  and  views  the  entire  United  States. 
The  infrared  (IR)  scanner  on  GOES  receives  outgoing  radiation  emitted  from 
the  earth's  surface  (and  any  emitting  mass  between  the  surface  and  the 
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Figure  8.     Required  inputs  for  the  P-model. 
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satellite),  then  calibrates  the  radiation  to  temperature  through  the 
Stephan-Boltzmann  Law.  The  best  pixel  resolution  for  GOES  is  3.6  km,  but 
due  to  distortion  from  the  earth's  curvature,  a  GOES  West  pixel  represents 
a  4.0  X  8.6  km  area  over  Florida. 

Good  agreement  between  satellite-sensed  temperatures  and  ground 
measurements  has  been  observed  for  clear  radiative  night  conditions. 
Martsolf  and  Gerber  (1981)  noted  ground  measurements  within  1  C  of 
satellite-sensed  temperatures.  Chen  et  al.  (1983)  compared  ground 
measurements  and  satellite  observations  from  clear  nights  during  winters 
of  1978  through  1981.  They  arrived  at  a  mean  correlation  coefficient  of 
0.87  and  an  average  standard  deviation  from  regression  of  1.57  C. 

BLM  Inputs 

The  BLM  required  a  sounding  of  the  atmospheric  boundary  layer  as 
input,  including  measurements  of  temperature,  wind  velocity  components, 
geostrophic  wind  components,  and  mixing  ratio  at  100-meter  increments  up 
to  3000  meters.    An  example  of  a  BLM  input  file  is  shown  in  Table  1. 

The  air  temperatures  were  converted  to  potential  temperatures  (Iribarne 
and  Godson,  1973): 


e 


T  (1000/P  )^, 


(15) 


where 


e 


potential  temperature  (K) 


T 


temperature  at  a  given  pressure  level  (K), 


atmospheric  pressure  (mb). 
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Table  1.  Sample  BLM  input  file, 
in  Table  2. 


Values  at  top  are  explained 


GAINESVILLE  SOUNDING 


TA=  10.00 
UGA  =  4.2U 
ZO  =  1.5000 
TX  =  40. 
TG  =  11.00 
RHOG  =  100.00 
SIGMAF  =  .75 


UA  =  4.24 
VGA  =  0.00 
CAPG  =  107193.560 
TN  =  0. 
RHOGX  =  100.00 
CSD  =  2928800. 


VA  =  0.00 
GLATD  =  29.5 
TM  =  11.00 
GOTIME  =  420.00 
RHOWLT  =  50.00 
CSW  =  4190000. 


QA  =  .0102 
DECLD  =-20.0 
IDOWN  =  319.0 
TRANSM  =  .95 
RHOM  =  75.00 
FN  =  3.0 


POTENTIAL  TEMPERATURES  T(I)  (C) 


8.858 

9.048 

9.105 

9.300 

9.527 

9.825 

10.044 

10.320 

11.205 

13.797 

14.473 

15.185 

15.957 

16.768 

18.342 

20.028 

21.636 

23.184 

24.621 

25.307 

26.091 

26.808 

27.655 

28.407 

29.175 

30.051 

30.887 

WIND  COMPONENTS 

U(I)  ( 

-K 

m  s  ) 

3.821 

4.702 

5.270 

5.838 

6.407 

6.975 

7.543 

8.112 

8.680 

9.816 

10.385 

10.953 

11.521 

12.117 

12.859 

13.598 

14.337 

15.076 

16.553 

17.292 

18.031 

18.769 

19.508 

20.247 

20.986 

21.724 

22.463 

WIND  COMPONENTS  V(I)      (m  s"b 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

MIXING 

RATIOS  Q(I)X1000 

(kg  water  vapor 

per  kg 

air) 

2.108 

2.003 

1.978 

1.640 

1.339 

1.133 

1.028 

.999 

.843 

.589 

.584 

.580 

.575 

.602 

.997 

1.671 

2.738 

4.140 

3.332 

2.982 

2.668 

2.385 

2.131 

1.903 

1.699 

1.515 

1.351 

GEOSTROPHIC  COMPS.  UG(I) 

(m  s" 

3.821 

4.702 

5.270 

5.838 

6.407 

6.975 

7.543 

8.112 

8.680 

9.816 

10.385 

10.953 

11.521 

12.117 

12.859 

13.598 

14.337 

15.076 

16.553 

17.292 

18.031 

18.769 

19.508 

20.247 

20.986 

21 .724 

22.463 

GEOSTROPHIC  COMPS.  VG(I) 

(m  s" 

') 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 
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X  =  .286,  or  R/Cp, 

R  =  ideal  gas  constant  (J  kg"^  K"^)» 

-2  -1 

C    =  specific  heat  (kJ  kg     C  ). 


The  potential  temperature  is  the  temperature  of  a  parcel  of  air  if  its 
ambient  volume  was  compressed  or  expanded  to  a  pressure  of  1000  mb  by  a 
dry  adiabatic  process.  This  is  done  as  a  normalization  of  the  temperature. 
One  reason  for  this  was  to  make  testing  of  the  model  easier.  For  example, 
the  potential  temperature  can  be  set  to  a  single  value  for  a  neutrally 
stable  boundary  layer  profile. 

The  mixing  ratio  is  the  fraction  of  mass  water  vapor  per  unit  mass 
dry  air.  It  was  calculated  from  the  dry  bulb  and  dew  point  temperatures 
using  the  following  relationships  (Iribarne  and  Godson,  1973;  Merva,  1975): 

log  U  =  -.000425  T  T^/(T-T^),  ^^^^ 
In  e    =  -2937. 4/T  -4.9283  log  T  +23.5518,  (17) 

3 

e  =  e  U,  (18) 

3  ' 

w  =  0.522  e/p,  (19) 
where 

T  =  dry  bulb  temperature  (K), 
T^  =  dew  point  temperature  (K), 
U  =  relative  humidity, 

e  =  vapor  pressure  at  ambient  temperature  (kpa), 

=  saturation  vapor  pressure  (kpa), 
w  =  mixing  ratio  (kg  water  vapor/kg  dry  air), 
p  =  atmospheric  pressure  (kpa). 
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The  geostrophic  wind  is  the  flow  that  would  be  produced  solely  by  the 
balance  between  the  pressure  gradient  force  and  the  coriolis  effect,  in 
the  absence  of  friction  (Sutton,  1953): 

\   =2F^7In*^HP  ' 

and , 

V„p   =     |-  +   |-      ,  (21) 
H  3  X        3  y 

where 

0)  =  angular  velocity  of  earth's  rotation  (s"^), 

4  =  geographic  latitude  (degrees  north), 

p  =  air  density  (kg  m~^), 

2w3in(l)  =  coriolis  parameter  (s~^), 

p  =  air  pressure  (n  m~^), 

t  =  time  (s) , 

x,y  =  horizontal  distances  (m) . 

The  measured  wind  near  the  surface  deviates  from  the  the  geostrophic 
flow  due  the  surface  frictional  effects.  The  measured  wind  approximates 
the  geostrophic  wind  within  and  above  the  higher  levels  of  the  boundary 
layer.    Mcintosh  (1972)  places  the  transition  level  at  about  600  meters. 

The  geostrophic  flow  was  calculated  from  the  measured  wind  by  use  of 
the  pressure  gradient  force  (Clarke  et  al.,  1971).  The  pressure  gradient 
was  determined  from  the  sounding  measurements  through  a  numerical 
approximation  of  the  Laplacian; 


Pl        P2   *   P3   -   P4   -  4p^ 

7|P    =   , 

b 


(22) 
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where 

p  =  pressure  of  surrounding  locations  (mb), 
p^  =  pressure  at  central  location  (mb), 

b  =  distance  from  central  station  to  surrounding  stations  (ra) . 

There  were  additional  inputs  to  the  3LM.  Several  of  these  depended 
on  location,  time  of  year,  and  soil  conditions.  Appropriate  values  for 
Florida  locations  were  used  (Table  1),  and  are  explained  in  Table  2.  .Some 
assumptions  were  made  in  lieu  of  actual  measurements.  The  high  atmospheric 
transraissivity  is  for  a  clear  dry  night.  The  field  capacity  and  wilt 
limit  are  for  sandy  soil  found  in  Florida  agricultural  regions.  The  soil 
moisture  and  slab  soil  moisture  are  assumed;  the  effects  of  estimation 
errors  are  shown  in  the  sensitivity  analysis  section. 

Procedure 

The  P-model  and  the  BLM  each  required  input  data  from  different 
outside  sources.  This  section  describes  the  input  data  sources,  acquisition 
of  the  data,  and  the  process  to  produce  the  temperature  predictions. 

Ground  Measurements 

Model  output  was  verified  with  1.5  m  air  temperature  measurements. 
These  were  recorded  at  Gainesville  and  Ruskin  by  different  means.  A 
station  similar  to  the  automated  weather  station  shown  in  Figure  8  was 
set  up  for  the  winter  of  1984  at  the  University  of  Florida's  Horticultural 
Unit.     Measurements  were  logged  on  a  Campbell  Scientific  CR5  Digital 
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Table  2.    Additional  BLM  inputs. 


Parameter  Value  Source 


Geographical  Latitude 

29.5  (degrees  North) 

Solar  Declination 

—ei\j,\j  (.aegrees; 

roughness  parameter 

1.5    (m)  (grove  environment) 

1 

vegetated  fraction  of  surface 

0.75  (grove  environment) 

assumed 

ground  albedo 

0.18 

2 

vegetative  albedo 

0.2 

2 

atmospheric  transmissivlty 

0.95  (clear  night  conditions) 
2928800  J  ra"^  C"^ 

assumed 

heat  capacity  (dry  soil) 

3 

soil  moisture — field  capacity 

100.0  (kg  m'^) 
50.0  (kg  m'^) 
75.0  (kg  m"^)  H^O 
100.0  (kg  m"^)  H^O 

4 

wilt  limit 

4 

subsoil 

assumed 

slab  soil  moisture 

assumed 

leaf  area  index 

3.0 

assumed 

Rosenberg  (1974) 
^Van  Wijk  (1963) 
^Armson  (1977) 
^Brady  (1974) 
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Recorder.  Temperatures  were  measured  by  copper-oonstantan  thermooouplea. 
The  thermocouples  were  calibrated  in  a  well-stirred  Ice  bath  and  the  bath 
temperature  was  recorded  over  several  minutes  by  the  CR5.  The  thermocouples 
averaged  0.6+0.1  C  over  the  time  period.  The  thermocouples  were  also 
compared  to  a  glass  thermometer  over  a  temperature  range  of  -8.0  to  +8.0 
C  in  a  well  stirred  salt  solution  (Figure  9).  The  P-model  had  been  running 
operationally  without  net  radiation  measurements  from  the  automated  weather 
stations  over  the  last  four  years,  so  no  net  radiometer  was  set  up  for 
this  study.  Measurements  were  made  in  the  field  at  the  levels  indicated 
in  Figure  8.    Measurements  were  taken  at  10-minute  intervals. 

Attempts  were  made  to  have  a  similar  station  operating  at  the  Weather 
Service  Office  (WSO)  in  Ruskin,  Florida,  during  the  Winter  of  1984-85, 
but  the  station  was  never  upgraded  to  an  operational  mode.  A  thermograph 
located  at  the  Ruskin  WSO  on  a  grass  surface  about  30  meters  from  any 
man-made  structure,  placed  1.5  meters  above  the  ground,  was  used.  The 
meteorologists  at  the  Ruskin  office  suggest  that  this  area  was  representative 
of  an  agricultioral  environment.  The  thermograph  was  checked  for  calibration 
weekly,  and  had  an  associated  measurement  error  of  +  1.0  C.  Temperature 
values  were  read  from  the  thermographs  at  30-minute  intervals. 

The  quick  response  time  of  thermocouples  and  the  greater  sampling 
frequency  of  the  CR5  provided  much  more  detailed  information  about  the 
temperature  changes  at  the  Gainesville  site  than  at  the  Ruskin  site.  The 
thermograph  tends  to  smooth  out  the  fluctuations,  as  shown  in  cooling 
curves  found  in  appendices  3  and  C. 
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Figure  9.     Calibration  of  CR5  with  thermocouples. 
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Satellite 

Satellite  temperatures  were  used  as  initial  surface  temperature  for 
the  modified  P-model.  Temperatures  of  nine  pixels  were  averaged  to 
determine  the  surface  temperature.  This  was  done  for  both  the  Gainesville 
and  the  Ruskin  sites. 

The  thermal  IR  data  were  acquired  by  the  Gainesville  HP  every  half 
hour  (Martsolf,  I982).  The  satellite  scans  the  entire  hemisphere  and 
transmits  the  raw  binary  data  to  Wallop's  Island,  Virgiaia,  where  gridding 
is  added.  The  stretched  data  are  retransmitted,  and  the  binary  data  are 
received  through  the  Gainesville  antenna.  The  data  are  reduced  into  a 
primary  sector  by  a  Schlumberger  EMR  822  sector izor,  and  a  Florida  sector 
consisting  of  an  array  of  I8O  x  100  pixels  is  produced  through  software. 

Slight  drift  in  the  satellite  caused  misalignment  between  maps.  The 
program  that  averages  the  nine  pixels  may  not  be  using  the  same  pixels 
from  one  map  to  the  next  if  the  maps  are  not  registered.  Therefore,  two 
alignment  routines  were  developed  to  search  fo"^  an  outstanding  map  feature. 
The  first  routine  used  the  map  grid  bits  to  find  the  outline  of  Lake 
Okeechobee.  The  grid  outline  of  Lake  Okeechobee  was  compared  to  the  grid 
outline  from  a  base  map.  The  difference  in  east-west  and  north-south 
pixel  locations  between  the  base  map  and  the  map  under  consideration 
determined  the  offsets.  A  second  routine  was  witten  that  searches  the 
actual  IR  data  for  Lake  Okeechobee.  The  lake  temperatures  are  usually 
warmer  than  the  surrounding  land  areas  during  a  frost  night.  The  program 
searched  for  and  highlighted  the  strongest  temperature  gradients,  producing 
an  outline  of  the  lake.  The  outline  was  then  compared  to  a  base  map 
outline  and  the  best  fit  was  calculated.    East-west  and  north-south  pixel 


44 


offsets  were  determined  and  added  to  the  latest  map  to  align  it  with 
the    first  map  of  the  night.    The  first  method  was  satisfactory  only  if 
the  gridding  was  placed  correctly  onto  the  map,  but  was  necessary  in  case 
there  were  clouds  over  Lake  Okeechobee  early  in  the  evening. 

The  average  of  nine  pixels  surrounding  each  automated  weather  station 
location  was  calculated  to  provide  a  substitiate  for  the  surface  temperature. 
Atmospheric  corrections  were  not  necessary  because  the  atmosphere  absorbs 
very  little  outgoing  radiation  on  clear  winter  nights  (Sutherland  et  al., 
1979).  Sach  pixel  provides  the  average  surface  temperature  of  the  34.0 
square  kilometer  area,  and  nine  were  averaged  to  smooth  anomalous  pixel 
temperatures . 

The  satellite  provided  surface  temperature  and  a  value  of  net  radiation 
^^ch  was  calculated  from  the  surface  temperature  through  the  Stephan-Boltzmann 
law.  Since  the  satellite  measured  only  outgoing  radiation,  the  calculated 
surface  term  ignored  the  back  radiation  from  the  sky,  and  therefore 
overestimated  surface  cooling  for  the  rest  of  the  night.  Surface  emissivity 
is  usually  in  the  range  of  0.90  to  0.98  (Sellers,  1965),  and  the  effective 
clear  sky  emissivity  is  between  0.5  and  0.7  (Monteith,  1980).  The  sky 
emissivity  is  directly  related  to  the  water  vapor  content  to  the  lower 
atmosphere,  and  increases  with  increasing  temperature  and  water  vapor. 
An  emissivity  of  0.66  was  used  in  the  equation  to  account  for  the  effects 
of  back  radiation. 

Atmospheric  Soundings 


The  National  Weather  Service  Automated  Field  Operational  System  (AFOS) 
provided  input  to  the  BLM.    Upper  air  soundings  from  Centreville,  Alabama, 
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Waycross,  Georgia,  Apalachicola,  Ruskin,  Cape  Canaveral,  West  Palm  Beach, 
and  Key  West,  Florida  provided  atmospheric  measurements  through  auid  beyond 
the  boundary  layer  at  00:00  and  12:00  CUT  daily. 

The  TTAA  atmospheric  sounding  report  (Table  3)  provided  dry  bulb  and 
dew  point  temperature,  height  of  the  pressure  surface,  wind  speed,  and 
wind  direction  at  mandatory  pressure  levels.  The  mandatory  levels  were 
surface,  1000  mb,  850  mb,  TOO  mb,  500  mb,  400  mb,  and  200  mb.  The  TTBB 
sounding  report  provided  pressure,  dry  bulb  temperature,  and  wet  bulb 
temperature  at  significant  levels  (Table  3). 

The  TTAA  sounding  reported  at  least  three  measurements  that  are  within 
the  boundary  layer:  surface,  1000  mb,  and  850  mb  levels.  The  700  mb 
level  is  usually  located  around  3000  meters.  Just  beyond  the  boundary 
layer.  Several  measurements  from  the  TTBB  report  often  fell  within  the 
first  3000  meters  of  the  atmosphere,  so  it  provided  additional  inputs. 
The  100-meter  interval  inputs  required  for  the  boundary  layer  model  could 
be  interpolated  from  these  sounding  measurements,  since  measuranents  from 
the  first  3000  meters  were  used.  A  linear  regression  was  first  run  on 
height  versus  pressure.  This  was  necessary  to  determine  the  pressure  at 
each  100-meter  interval  so  that  each  potential  temperature  could  be 
calculated.  The  correlation  coefficients  of  10  pressure-vs. -height 
regressions  from  soundings  were  averaged.  The  results  gave  an  average  r 
of  .999  +  .001,  indicating  strong  linearity  between  the  pressure  and  height 
(Little  and  Hills,  1978). 

The  dry  bulb  temperature,  dew  point  temperature,  and  wind  speed  were 
extracted  from  the  sounding  through  the  use  of  piecewise  linear  interpolation 
(Burden  et  al.,  1978).  The  TTBB  report  indicates  changes  of  linear-ity  in 
dry  bulb  or  dew  point  vs.  height,  so  this  interpolating  procedure  should 
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Table  3.  Example  of  an  upper  air  sounding  from  AFOS. 


MIAMANTBW 
WOUSOO  KTBW  101200 

72210  TTAA  50121  72210  99017  22205  09002  00161  25243  16504 
85579  13256  30503  70211  06846  06007  50591  07980  08012  40760 
20980  04508  30957  35359  34512  25091  455//  33510  20235  551// 
34520  15416  651//  04046  10660  677//  06521  88124  703//  05536 
77999  51515  10164  00001  10194  26504  03504= 

MIASGLTBW 
WOUSOO  KTBW  161200 

72210  TTBB    5612/  72210  00017  22012  11000  24413  22904  19250 

33850  15838  44756  10460  55700  07450  55690  07058  77679  05840 

88653  04680  99518  06163  11500  08944  22494  09902  33453  12922 

44433  15155  55413  16350  66410  17560  77400  18756  88388  19759 

99375  21558  11363  22956  22300  33580  33274  38964  44190  589// 
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provide  good  representation  of  the  actual  interval  values  (Cahir  et  al., 
1978).  The  complete  input  for  the  BLM  was  produced  through  the  following 
sequence: 

a.  Read  TTAA,  TTBB  reports,  collect  pressure  versus  height,  dry  bulb 
and  dew  point  temperatures,  and  wind  speed; 

b.  Sort  the  TTAA  and  TTBB  results  in  height-increasing  order; 

c.  Perform  linear  regression  on  pressure  vs.  height,  equation  calculated; 

d.  Perform  piecewise  linear  interpolation  for  temperature  vs.  height, 
temperatures  interpolated  at  100-m  intervals  starting  at  100  m; 

e.  Perform  piecewise  linear  interpolation  for  dew  point  temperature; 

f.  Perform  piecewise  linear  interpolation  for  wind  speed; 

g.  Calculate  potential  temperature  from  c  and  d; 

h.  Calculate  mixing  ratio  from  e; 

i.  Calculate  geostrophic  wind  speed  from  pressure  gradients. 

A  diagram  of  source  data  acquisition  and  transfer  is  shown  in  Figure 
10.  The  sounding  information  is  transferred  from  the  AFOS  Data  General 
NOVA  computer  to  an  HP  21MX-M  computer  at  the  Ruskin  National  Weather 
Service  Office.  The  data  are  moved  from  Ruskin  to  the  SFFS  HP  21MX-E  via 
a  DS-1000  link  over  a  dedicated  phone  line.  The  sounding  data  are  processed 
into  BLM  input  files  on  the  SFFS  computer.  The  input  files  are  then  sent 
to  the  Florida  Agricultural  Service  and  Technology  CYBER  730.  The  P-model 
was  run  on  the  HP  21MX-E,  the  BLM  on  the  CYBER  730.  Both  were  run 
interactively  from  user  sessions.  Details  of  the  procedure  to  run  the 
models  and  programs  used  to  process  the  input  files  are  found  in  Appendix 
E. 
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Figure  10.     Diagram  of  the  acquisition,  transfer,  and  processing 
of  data  for  model  predictions. 


RESULTS 


The  nights  observed  for  this  study  are  grouped  into  three  categories. 
The  first  is  the  radiative  nights,  during  which  the  dominant  cooling 
mechanism  in  the  surface  air  layer  was  radiative  heat  loss  to  the  sky. 
The  second  is  advective  nights,  when  cold  air  flow  into  the  area  was  the 
dominant  cooling  mechanism.  The  third  consists  of  episodes  that  didn't 
fit  into  either  catagory.  These  episodes  involved  situations  where  the 
microclimate  was  influenced  by  larger  scale  atmospheric  occurrences. 

Two  sites  are  analyzed  to  test  the  BLM  and  P-model  predictions  against 
observed  temperature.  Eleven  nights  are  observed  for  Ruskin  (RUS)  and 
twelve  nights  are  observed  for  Gainesville  (GNV)  (Table  4).  The  Boundary 
Layer  Model  and  Modified  P-model  were  run  on  both  the  Gainesville  and 
Ruskin  sites.  The  original  P-model  was  run  on  the  Gainesville  site  only, 
since  the  measurement  tower  was  not  put  into  operation  by  the  WSO  staff 
at  Ruskin  during  the  Winter  of  1984-85. 

Linear  regression  was  run  on  predicted  versus  observed  temperatures. 
Linear  regression  is  a  standard  approach  to  evaluating  model  performance, 
and  can  show  details  of  the  performance,  with  limitations.  Ideally,  the 
predicted  versus  observed  curve  should  have  an  r  value  of  1.0,  indicating 
a  straight  line,  and  a  slope  of  1.0,  indicating  the  difference  between 
the  observed  and  predicted  temperatures  did  not  change  at  each  time  period. 

Willmott  (1984)  reports  that  the  correlation  coefficient  r  and 
explained  vsu?icince  r   alone  may  not  be  indicative  of  actual  model  performance 
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Table  4.    List  of  nights  observed  for  this  study. 


site  radiative  advective  anomalous 


GNV 

Jan. 

6, 

1985 

Jan. 

5, 

1985 

Jan. 

13, 

1985 

Jan. 

16, 

1985 

Jan. 

21, 

1985 

Jan. 

22, 

1985 

Jan. 

20, 

1985 

Jan. 

26, 

1985 

Jan. 

23, 

1985 

Jan. 

24, 

1985 

Jan. 

27, 

1985 

Feb. 

10, 

1985 

Feb. 

17, 

1985 

BUS 

Dec. 

8, 

1984 

Dec. 

7, 

1984 

Jan. 

23, 

1985 

Dec. 

9, 

1984 

Jan. 

21, 

1985 

Jan. 

6, 

1985 

Jan. 

26, 

1985 

Jan. 

16, 

1985 

Jan. 

24, 

1985 

Jan. 

27, 

1985 

Feb. 

16, 

1985 
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under  certain  circumstanoes .  For  example,  the  predicted  versus  observed 
curve  may  give  an  r  value  of  1.0,  indicating  a  straight  line,  but  the 
slope  may  deviate  considerably  from  1.0  and  the  intercept  may  be  different 
from  zero.  In  this  case  the  model  performance  could  be  quite  poor. 
Willmott  suggested  other  methods  of  evaluating  model  performance.  Root 
mean  square  error  (RMSE)  and  Mean  Absolute  Error  (MAE)  are  two  comparison 
techniques  used  in  this  project. 

The  RMSE  and  MAE  represent  average  magnitude  differences  between  the 
observed  and  predicted  values.  This  is  a  very  useful  method  of  presenting 
the  error  associated  with  the  predictions.  The  RMSE  tends  to  be  greater 
than  the  MAE  because  the  differences  are  squared,  so  any  large  errors  will 
amplify  the  average  difference. 

T-tests  were  run  on  the  difference  between  predicted  and  observed 

values  (AT  =  T    -  T  )  to  determine  if  AT  is  significantly  different  from 
P  0 

zero  at  each  hourly  time  interval.  The  t-tests  were  run  at  a  95%  confidence 
level. 

Model  results  are  summarized  in  Table  5.  The  graphs  of  predicted 
versus  observed  temperatures  for  the  BLM,  P-model,  and  Modified  P-model 
for  all  radiative  nights  at  the  Gainesville  site  are  shown  in  Figures  11, 
12,  and  13>  respectively.  The  results  of  the  t-tests  for  the  BLM,  P-model, 
and  Modified  P-model  are  shown  in  Figures  14,  15,  and  16,  respectively. 
The  graphs  of  predicted  versus  observed  temperatures  and  results  of  t-tests 
for  all  radiative  nights  at  the  Ruskin  site  are  shown  in  Figures  17-19. 
The  predicted  versus  observed  temperatures  for  the  BLM  and  P-model  during 
advective  nights  for  the  Gainesville  and  Ruskin  sites  are  shown  in  Figures 
20-22. 
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Table  5.  Summary  of  model  results  from  radiative  nights.    The  values 


shown  are  the  average  and  standard  deviations  of  95%  confidence 
intervals  from  t-tests  run  on  the  hypothesis  that  the 
difference  between  predicted  and  observed  temperatures  is 
not  significantly  different  from  zero.  The  averages  are 
for  all  hours  of  radiation  and  advection  nights  from  each 
site. 


site 


model 


95%  confidence  interval  (C) 


GNV 


Mod.  P-model 


BLM 


P-model 


3.4  +  1.2 
5.8  +  1.0 
5.7  +  0.9 


RUS 


Mod.  P-model 


BLM 


2.7  +  0.6 
18.8  +  2.U 


•*Mod.  P-model 


9.2  +  1.5 


**cloudy  nights  removed 
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Figure  11. 


BLM  predicted  versus  observed  temperatures  for 
all  radiation  nights,  Gainesville  location. 
Solid  line  indicates  perfect  fit. 
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Figure  12, 


P-model  predicted  versus  observed  temperatures 
for  all  radiation  nights,  Gainesville  location. 
Solid  line  indicates  perfect  fit. 
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Figure  13.     Modified  P-model  predicted  versus  observed 

temperatures  for  all  radiation  nights,  Gaines- 
ville location.  Solid  line  indicates  perfect 
fit. 
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Figure  14.     Results  of  t-tests  on  difference  between 
BLM  predicted  and  observed  temperatures 
being  significantly  different  from  zero, 
showing  95%  confidence  intervals,  Gaines- 
ville location. 
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Figure  15.     Results  of  t-tests  on  difference  between 

P-model  predicted  and  observed  temperatures 
being  significantly  different  from  zero, 
showing  95%  confidence  intervals,  Gaines- 
ville location. 
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Figure  16.     Results  of  t-tests  on  difference  between 
Modified  P-model  predicted  and  observed 
temperatures  being  significantly  differ- 
ent from  zero,  showing  95%  confidence  in- 
tervals, Gainesville  location. 
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Figure  17.    BLM  predicted  versus  observed  temperatures  for 
all  radiation  nights,  Ruskin  location.  Solid 
line  indicates  perfect  fit. 
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Figure  18.     Modified  P-model  predicted  versus  observed  temperatures 
for  all  radiation  nights,  Ruskin  location. 
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Figure  19.     Results  of  t-tests  on  difference  between 
BLM  predicted  and  observed  temperatures 
being  significantly  different  from  zero, 
showing  95%  confidence  intervals,  Ruskin 
location. 
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Figure  20.     BLM  predicted  versus  observed  temperatures  for 
all  advection  nights,  Gainesville  location. 
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Figure  21 o    BLM  predicted  versus  observed  temperatures  for 
all  advection  nights,  Ruskin  location. 
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Figure  22.    P-model  predicted  versus  observed  temperatures  for 
all  advection  nights,  Gainesville  location. 


DISCUSSION 

The  observed  cooling  curve  followed  a  decaying  exponential  curve, 
typical  for  a  radiative  frost  night,  in  six  of  the  twelve  nights  studied 
at  the  Gainesville  location  and  seven  of  the  eleven  nights  studied  at  the 
Ruskin  site.  Of  the  remaining  nights,  three  of  the  Gainesville  nights 
and  three  Ruskin  nights  were  advective,  and  three  Gainesville  nights  and 
one  Ruskin  night  showed  unusual  meteorological  circumstances  that  fall 
into  an  anomalous  category.  Discussion  of  the  individual  nights  and 
cooling  curves  from  the  individual  nights  are  in  Appendices  B  and  C. 

Radiative  Nights 

The  typical  radiative  night  cooling  curve  is  produced  from  strong 
radiative  heat  loss  from  the  ground  and  vegetative  surfaces,  and  from  the 
air  in  the  layer  near  the  ground  surface.  A  divergence  of  heat  occurs  at 
the  ground  surface  when  the  radiative  heat  loss  is  greater  than  the 
conduction  of  heat  through  the  soil  to  the  soil  surface.  The  radiative 
heat  flux  decreases  as  the  surface  cools.  Blackadar  (  1979)  suggested 
another  contributing  factor  to  the  decaying  exponential  curve.  The  wind 
shear  between  the  upper  layers  and  the  surface  air  layer  is  small  in  the 
early  evening  because  the  upper  level  winds  have  Just  begun  to  accelerate. 
As  the  upper  winds  increase,  so  does  the  wind  shear,  which  in  turn  increases 
mixing  at  the  lower  levels. 
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Linear  regressions  were  run  starting  with  the  fourth  forecast  hour 
when  comparing  the  P-model  to  the  BLM.  The  reason  for  this  is  that  the 
first  three  hours  of  the  observed  values  form  the  base  for  the  P-model 
predictions,  and  therefore  the  P-model  curve  is  the  same  as  the  observed 
for  those  hours.  The  predictions  then  begin  to  deviate  from  the  observed 
after  the  third  hour. 

Gainesville 

The  y  intercept  of  the  linear  regression  curve  indicates  that  the 

BLM  predicted  the  cooling  curve  warmer  than  the  observed  (Figure  11,  page 

53).    This  is  discussed  in  the  sensitivity  analysis  section.  However, 

based  on  the  results  of  the  t-tests,  the  variance  between  predicted  and 

observed  temperatures  was  considerably  less  for  the  BLM  than  for  the 

P-model  throughout  the  night.    The  BLM  curves  were  very  similar  in  fit  to 

the  observed  curves  under  ordinary  radiative  conditions  as  shown  by  a 

2 

slope  of  1.09  and  r  of  .72,  although  displaced  upwards  by  an  average  of 
3«35  C.  The  difference  between  the  observed  and  predicted  is  significantly 
different  from  0.0  at  most  hours  for  the  BLM  (Figure  14,  page  56),  because 
of  the  3.35  C  displacement.  The  95%  confidence  interval  resulting  from 
the  t-tests  averages  +  1.7  C  over  the  night.  The  difference  between 
predicted  and  observed  temperatures  for  the  P-model  varies  much  more, 
predicting  too  warm  in  some  cases,  and  too  cold  in  others.  The  slope  of 
the  regressed  temperatures  is  0.88,  (Figure  12,  page  54),  but  the  graph 
indicates  a  large  variation  in  difference  between  the  predicted  and  observed 
temperatures.  The  95%  confidence  intervals  average  +  2.9  C  over  the  night 
for  the  P-model  (Figure  15,  page  57).    The  Modified  P-model  results  were 
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similar  to  the  P-model  results  (Figure  13,    page  55)  with  confidence 
intervals  averaging  +  2.9  C  (Figure  15,  page  58). 

Ruskin 

The  BLM-predicted  cooling  waa  very  close  to  the  observed  cooling  for 

radiation  nights  at  the  Ruskin  sounding  site.     The  slope  was  0.9,  the 

2 

intercept  value  -.2  C,  and  r  was  0.7  (Figure  17,  page  59).  The  modified 
P-model  produced  large  errors  in  its  predictions  at  the  Ruskin  site,  due 
to  the  presence  of  clouds  during  the  prediction  base  hours  (Figure  18, 
page  60).  The  satellite  receives  radiation  from  the  cold  cloud  tops,  and 
as  the  clouds  clear  the  satellite-measured  temperatures  begin  to  increase. 
The  increase  in  base  temperature  forced  the  P-model  to  predict  warming 
instead  of  cooling.  This  occurred  on  two  of  the  seven  radiative  nights 
at  Ruskin.  Excluding  those  two  nights,  the  modified  P-model  still  had 
large  variation  in  difference  between  predicted  and  observed  temperatures 
(95$  confidence  intervals  of  +  4.6  C).  The  t-tests  showed  that  the 
BLM-predicted  temperatures  for  the  Ruskin  site  were  closer  to  the  observed 
temperatures  than  for  any  other  case  (Figure  19,  page  51).  There  was  no 
significant  difference  between  the  predicted  minus  observed  temperature 
difference  and  zero  for  11  of  the  13  time  periods,  and  the  95%  confidence 
interval  value  averaged  +  1.4  C. 

These  results  show  that  the  BLM  predicted  more  precisely  than  the 
P-model  for  the  radiative  nights  at  both  sites.  This  is  due  to  the  BLM 
modeling  three  important  processes  in  more  depth  than  the  P-model:  The 
surface  radiation  flux,  the  turbulent  heat  transfer,  and  presence  of 
vegetation.    First,  the  initial  radiation  terms  of  the  P-model  are  either 
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assumed  or  calculated  from  the  measured  ground  surface  temperature.  An 
exponent:.al  regression  is  then  performed  using  the  three-hour  base  radiation 
measurements.  The  assumption  that  the  radiation  flux  will  decay  in  an 
exponential  manner  is  good,  but  using  assumed  values  or  calculating 
radiation  flux  from  surface  temperatures  (ivhich  are  inherently  difficult 
to  measure  accurately  with  thermocouples),  will  tend  to  lessen  the  accuracy 
of  the  radiation  values.  The  BLM,  on  the  other  hand,  calculates  the 
radiation  term  from  a  ground  surface  temperature  that  has  been  determined 
from  calculations  of  soil  heat  flux  and  radiation  flux  divergence.  The 
new  radiation  terms  are  calculated  at  each  time  step  (two  minutes). 

Another  process  that  is  handled  better  by  the  BLM  is  the  turbulent 
convective  heat  transfer.  The  P-model  considers  only  the  first  9  m  of 
the  boundary  layer,  where  the  wind  profile  decays  exponentially  from  the 
9  m  velocity  to  zero  at  z^.  Neither  vertical  heat  nor  momentum  transfer 
is  accounted  for  by  the  P-model.  Only  a  convective  heat  transfer  coefficient 
is  calculated  (see  Models  and  Methods  chapter).  The  BLM  takes  into  account 
the  heat  and  momentum  transfer  from  the  air  layers  above  the  surface  air 
layer,  and  calculates  a  new  turbulent  eddy  transfer  coefficient  at  every 
time  increment. 

Finally,  the  BLM  includes  a  vegetation  model;  the  P-model  does  not. 
The  vegetative  canopy  acts  as  a  radiation  absorber  and  emitter.  The 
vegetation  section  also  accounts  for  ground  soil  and  plant  moisture  whi.ch 
can  add  to  the  evaporative  effects  as  well  as  instability  in  the  surface 
air  layer,  although  these  are  small  at  night. 
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Advective  Nights 

Cold  air  was  rapidly  moving  into  the  prediction  eirea  with  the  passage 
of  a  cold  front  on  certain  nights.  Since  the  soundings  used  for  BLM  input 
were  taken  at  a  single  point  in  time,  the  impending  arrival  of  the  advected 
cold  air  mass  was  not  always  taken  into  consideration  by  the  model.  The 
magnitude  of  the  BLM  predicted  and  observed  temperature  difference  was 
not  large  for  the  Gai  nesville  site  (Figure  20,  page  62),  but  the  model 
underpredicted  the  cooling.  The  difference  increased  as  the  nights 
progressed.  The  difference  for  Ruskin  was  much  greater  (Figure  21,  page 
63)  as  Indicated  by  a  y  intercept  of  9.2  C.  Most  of  this  error  was  from 
the  night  of  January  21. 

The  errors  in  the  P-model  predictions  were  slightly  smaller  than  the 
BLM  for  Gainesville  advective  nights.  The  temperature  trend  could  have 
been  detected  by  the  P-model  if  the  advection  was  occurring  during  the 
three-hour  base  period  (Figure  22,  page  64),  making  the  predictions  closer 
to  the  observed. 

Anomalous  Nights 

In  certain  situations  meso-  and  synoptic-scale  processes  were  occurring 
that  the  models  could  not  anticipate.  These  processes  tended  to  moderate 
the  cooling.  The  BLM  predicted  the  worst  case  given  an  initial  set  of 
conditions  in  these  situations,  so  the  moderation  prevented  the  actual 
temperatures  from  dropping  below  the  predicted  temperatures. 
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Gainesville  BLM  Predictions  Compared  to  Observed  Satellite  Data 

The  Gainesville  predictions  were  compared  to  the  satellite  temperatures 
since  the  satellite  covers  a  larger  area  and  therefore  is  spatially 
averaged.  Satellite  temperatures  and  BLM  predictions  versus  time  are 
shown  in  Figures  23-29.  The  BLM  did  not  predict  the  satellite  tenperatures 
very  accurately.  The  only  night  that  the  BLM  predictions  were  similar  to 
the  satellite  ivas  on  the  2Uth  of  January  (Figure  27).  The  reason  for  this 
may  be  that  the  near-surface  temperature  can  vary  to  a  great  extent 
horizontally  during  a  frost  night,  depending  on  surface  features  and 
topography  (Bartholic  and  Martsolf ,  1979).  Chen  et  al.  (1983)  showed 
that  the  temperatures  across  Florida  counties  can  vary  by  several  degrees 
Celsius. 

Sensitivity  Tests 

Sensitivity  tests  were  performed  on  the  BLM  to  determine  the  cause 
of  the  Gainesville  predictions  being  too  high.  The  temperature  .predictions 
from  the  night  of  February  15-17  (0M8)  were  tested  for  sensitivity  to 
various  inputs. 

The  sensitivity  analysis  used  the  following  equations: 

absolute: 

s(y^ ,p) 

relative: 

S(y^ ,p) 
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Figure  23,     Average  and  standard  deviations  of  nine  satellite 
pixels  and  BLM  temperature  predictions  versus 
time,  Gainesville  location,  Jan.  6,  1985. 
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Figure  24.    Average  and  standard  deviations  of  nine  satellite 
pixels  and  BLM  temperature  predictions  versus 
time,  Gainesville  location,  Jan.  16,  1985. 
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Figure  25.    Average  and  standard  deviations  of  nine  satellite 
pixels  and  BLM  temperature  predictions  versus 
time,  Gainesville  location,  Jan.  22,  1985. 
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Figure  26. 


Average  and  standard  deviations  of  nine  satellite 
pixels  and  BLM  temperature  predictions  versus 
time,  Gainesville  location,  Jan.  23,  1985. 
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Figure  27. 


Average  and  standard  deviations  of  nine  satellite 
pixels  and  BLM  temperature  predictions  versus 
time,  Gainesville  location,  Jan.  24,  1985, 
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Figure  28. 


Average  and  standard  deviations  of  nine  satellite 
pixels  and  BLM  temperature  predictions  versus 
time,  Gainesville  location,  Jan.  27,  1985. 
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Figure  29,    Average  and  standard  deviations  of  nine  satellite 
pixels  and  BLM  temperature  predictions  versus 
time,  Gainesville  location,  Feb.  17,  1985. 
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where 

y  =  the  output  variable  to  be  analyzed, 

p  =  the  input  variable  or  parameter  to  be  analyzed. 

Equation  23  indicates  the  percentage  change  in  output  variable  y  given  a 
percentage  change  in  input  parameter  p.  The  surface  air  layer  temperature 
is  the  output  variable  studied  here.  The  percentage  of  change  in  the 
inputs  was  based  on  the  magnitude  of  error  that  was  associated  with  the 
input  measurements  or  variability  of  an  assumed  parameter.  For  example, 
the  error  associated  with  the  satellite  temperature  measurements  was  on 
the  order  of  +  1.0  C,  so  the  effect  of  a  ^%  change  in  the  input  surface 
temperature  was  used. 

The  results  of  sensitivity  tests  are  shown  in  Tables  6  and  7.  Two 
forecast  hours,  the  fourth  and  the  twelth,  a^e  selected  to  show  if  the 
effects  of  the  change  in  initial  input  is  increasing  or  decreasing  the 
change  in  output  as  the  night  progresses. 

The  results  show  that  the  surface  air  temperature  is  most  sensitive 
to  the  initial  soil  properties  and  temperatures,  and  the  initial  surface 
air  layer  temperature.  The  soil  properties  affect  the  soil  heat  flux 
which  in  turn  affects  the  soil  surface  temperature  T^.  The  outgoing 
radiation  is  determined  from  T^.  Since  radiation  is  the  dominant  process 
in  the  first  meter  of  the  surface  air  layer,  these  inputs  would  have  the 
greatest  effect  on  the  air  surface  temperature.  Later  in  the  night  the 
variables  affecting  the  turbulent  heat  exchange  become  important,  as  shown 
by  the  increase  in  sensitivity  to  and  z^,  but  the  soil  properties  still 
dominate.  Note  the  strong  effects  of  the  initial  temperatures;  only  a  ^% 
change  In  T  ,  T  ,  and  T    brings  about  a  .1  to  .6%  change  in  predicted  T  . 


Table  6.    Inputs  and  parameters  tested  in  sensitivity  analysis. 


paraine  ter 

0  vmhol 

innut  value 

sfc  air  temp 

\LaJ 

280.9  K 
U  77  m 

^•11     ui  0 

U  77  m  «^"^ 

sfc  air  wind  speed 

f  Ma  ■) 

sfc  air  geostrophic  wind  speed 

\.  uga ; 

roughness  length 

1    0  m 

in71QU  0  J  m"^ 

slab  heat  capacity 

(.capg; 

base  Temp 

(Tm) 

281 .9  K 

atmospneric  transmisaxviuy 

(  V  r*fl  n  ^Tn  ^ 

ground  sfc.  temp 

(Ts) 

281 .9  K 

soil  iieiu  capacity 

100.0  ke  m 

3011  Wilt  iimib 

50.0  kg  m"^ 

3UD3011  Wauer  uonteno 

75.0  kg  m"^ 

( fhOff ) 

50.0  kg  m"^ 

(csd) 

2.0x10e+6  J  m"^  C~ 

water  heat  capacity 

(csw) 

H.  lyxiue+Q  0  UI  L» 

leaf  area  index 

(LAI) 

3.0 

vegetative  coverage 

(sigmaf ) 

0.75 

sfc  air  mixing  ratio 

(qa) 

.0076  kg/kg 

daily  IR  cooling 

(rcool) 

. 833x1 Oe-3  C  hour" 
1 . 0  cm  hour" 

subsidence 

(wdown) 

mixing  depth 

(Zl) 

60.0  m 

.66614  W  m"^  C"^ 
.5  ms 

soil  thermal  conductivity 

(lamdag) 

turbulent  heat  trans,  coef. 

(K) 
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Table  7.    Results  of  sensitivity  tests.    Ap  is  the  change  in  the 
input  parameter  given  in  equations  22  and  23. 


relative  sensitivity  Ta  at  hour 


%  change 
input 

iJchange  Ta 
4.00  12.00 

4.00 

+Ap 

-Ap 

12. 
+Ap 

00 

-Ap 

Ta 

1.0 

.173 

.066 

4  18 

-.31 

-.49 

Ua 

100.0 

0.00 

0.00 

ii  UP 

4.42 

-.40 

-.40 

Uga 

100.0 

0.00 

0.00 

U  UP 

4.42 

-.40 

-.40 

zo 

20.0 

.079 

.088 

u  81 

5-03 

-.01 

.23 

capg 

10.0 

0.00 

0.00 

4.42 

-.40 

-.40 

Tm 

1 .0 

.299 

.641 

4  82 

^.Q9 

.47 

-1 .28 

transm 

10.0 

0.00 

0.00 

4  42 

4.42 

-.40 

-.40 

Tg 

1.0 

.270 

.095 

ll  7Q 

4  04 

-.27 

-.53 

rhogx 

20.0 

.054 

.264 

4  45 

4.33 

-.23 

-.59 

rhowlt 

20.0 

0.00 

0.00 

U  4? 

4.42 

-.40 

-.40 

rhom 

20.0 

0.00 

0.00 

4  42 

4.42 

-.40 

-.40 

rhog 

100.0 

.170 

.312 

7.  QQ 

-.54 

csd 

100.0 

.259 

.799 

4.70 

3.98 

.42 

-1.76 

csw 

10.0 

.072 

.015 

4.43 

4.42 

-.39 

-.41 

lai 

100.0 

-.123 

-.132 

4.35 

4.69 

-.47 

-.11 

sigmaf 

20.0 

-.068 

-.070 

4.34 

4.53 

-.48 

-.29 

qa 

100.0 

.274 

.990 

4.79 

4.03 

-.27 

-.54 

Zl 

100.0 

.022 

.579 

4.27 

4.21 

.01 

-1.57 

lamdag 

100.0 

.316 

.981 

4.72 

3.93 

.48 

-1.93 

rcool** 

100.0 

-.162 

-.338 

3.97 

4.42 

-1.32  . 

-.40 

wdown** 

100.0 

0.00 

.022 

4.65 

4.65 

.14 

.08 

K 

100.0 

.040 

.253 

4.47 

4.36 

-.12 

-.81 

»*These  two  parameters  were  set  to  0.00  in  the  original  model  run. 
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The  geostrophic  wind  has  an  Indirect  effect  on  the  temperature 

predictions,  but  is  very  important  in  the  development  of  the  predicted 

wind  profile.    Since  f(V  -V)  is  an  additive  term  (equations  4  and  5),  the 

S 

presence  of  geostrophic  wind  increases  the  predicted  wind  velocity  in  the 

lower  levels,  causing  the  predicted  low-level  nocturnal  jet  to  develop. 

The  measured  wind  tends  to  become  more  geostrophic  with  height,  so  the 

f(V  -V)  term  becomes  smaller. 
S 

The  geostrophic  wind  must  be  calculated  from  good  pressure  gradient 
data.  The  addition  of  geostrophic  wind  to  the  BLM  infile  often  caused  a 
fluctuation  in  the  predicted  temperature  similar  In  appearance  to  the 
fluctuations  observed  in  the  actual  environment.  However,  the  variation 
in  the  predictions  was  due  to  programmatic  numerical  instability  rather 
than  modeled  meteorological  processes.  Apparently,  the  spatial  density 
of  the  NWS  sounding  station  network  was  not  sufficient  for  calculating 
accurate  pressure  gradients,  therefore  the  geostrophic  determinations  were 
insufficient. 

The  idea  of  obtaining  calculated  geostrophic  wind  from  the  synoptic 
models  was  suggested  near  the  end  of  the  data  collection  period.  The 
National  Weather  Service  passes  a  file  containing  results  of  LFM  analysis 
and  predictions  through  AFOS.  This  file  (FRH63)  includes  the  geostrophic 
wind  for  the  lowest  50  mb,  from  the  analysis  time  through  forecasts  up  to 
48  hours  in  6-hour  forecast  increments. 

The  sensitivity  tests  do  not  reveal  any  particular  input  variable 
being  responsible  for  the  high  Gainesville  predictions.  This  led  to  the 
conclusion  that  the  consistently  high  BLM  predictions  for  the  Gainesville 
site  were  due  to  the  interpolation  of  the  input  sounding.  The  Gainesville 
sounding  was  produced  from  Interpolation  between  the  Huskin  and  Waycross 
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sounding  sites.  The  original  plan  was  to  interpolate  between  these  two 
plus  two  additional  soundings,  Apalachicola  and  Cape  Canaveral.  The  Cape 
Canaveral  site  does  not  follow  the  standard  NWS  sounding  schedule  so 
soundings  were  usually  unavailable  at  0000  CUT;  this  eliminated  the  use 
of  four  soundings.  Spatially,  Gainesville  is  located  almost  midway  between 
Waycross  and  Ruskin.  However,  Waycross  and  Gainesville  are  inlaind  sites, 
whereas  Ruskin  is  located  within  20  km  of  the  Tampa  Bay.  The  interpolated 
input  temperature  may  be  biased  upward  toward  the  Ruskin  profile,  when  in 
fact  the  boundary  layer  temperature  profile  for  Gainesville  may  be  better 
represented  by  the  Waycross  sounding. 

One  of  the  biggest  drawbacks  of  using  the  BLM  on  rioAr-ai- i.onal  basis 
is  the  problem  of  inputs.  The  model  does  very  well  at  the  sounding  site, 
but  the  sounding  network  is  not  dense  enough  to  cover  inland  areas  of 
Florida.  Interpolation  does  not  seem  to  be  the  best  method  of  obtaining 
data  between  sounding  stations.  The  possibility  exists  that  the  VAS  (VISSR 
atmospheric  sounder)  may  provide  sufficient  information  from  satellites 
to  develop  BLM  inputs  from  satellite  data  (Jedlovec,  1985).  This  would 
enable  the  model  to  predict  for  any  location  within  the  satellite  view, 
without  the  need  to  interpolate  between  radiosonde  launch  sites. 

The  VAS  radiometer  views  12  channels  to  discern  vertical  moisture 
and  temperature  profiles  through  the  atmosphere.  The  spatial  resolution 
is  similar  to  the  GOES  IR  (7.0  km),  as  opposed  to  the  low  density  sounding 
station  network  (100-200  km  over  land  and  practically  absent  over  the 
ocean).  Also,  the  potential  exists  for  temporal  increases  in  the  profiles, 
perhaps  as  often  as  every  half  hour.  This  could  enable  the  BLM  to 
self -correct  as  time  progresses  into  the  night. 
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The  VAS  data  are  being  researched  to  determine  the  VAS  accuracy  with 
respect  to  the  radiosonde  data.  VAS  data  are  difficult  to  obtain  on  an 
operational  basis,  but  in  the  future  this  information  could  become  quite 
valuable  as  model  input. 


SUMMARY  AND  CONCLUSIONS 


Two  atmospheric  models  were  analyzed  and  compared  to  determine  which 
would  more  accurately  predict  the  cooling  of  a  vegetated  environment  under 
nocturnal  conditions.  The  two  models  were  the  P-model,  which  was  developed 
on  the  Satellite  Frost  Forecast  System,  and  a  Boundary  Lfv>»r»  Model  (BLM) 
developed  by  Alfred  Blackadar. 

Atmospheric  soundings  were  used  as  input  to  the  BLM.  These  were 
obtained  from  National  Weather  Service  AFOS,  through  a  DS-1000  link  between 
Gainesville  and  the  WSO  in  Ruskin,  Florida.  Satellite  data  and  automated 
weather  station  measurements  were  used  for  P-model  input.  The  satellite 
data  were  obtained  from  GOES  West  through  a  direct  digital  downlink. 

Model-produced  predictions  were  analyzed  for  12  nights  at  the 
Gainesville  site  and  11  at  the  Ruskin  site.  The  predicted  versus  observed 
temperatures  were  compared  using  linear  regression,  t-test,  RMSE,  and  MAE. 
A  sensitivity  analysis  was  performed  on  the  Gainesville  site  for  the  night 
of  the  February  16-17,  1985. 

The  Boundary  Layer  Model  predicted  the  cooling  curve  more  precisely 
than  the  P-model  and  Modified  P-model  during  radiative  nights  in  the 
vegetated  environment.  For  Gainesville,  the  overall  regressions  of  observed 
versus  predicted  temperatures  from  the  fourth  forecast  hour  through  the 
remainder  of  the  night  were: 

3LM  predicted  =  1 .09*observed  +3.35      r    =  .72, 

P-model  predicted  =  0.88»observed  -  1.18      r^  =  .28, 

modified  P-model      predicted  =  0.69»ob3erved  +  0.93      r^  =  .27. 
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The  Ruskin  regressions  were: 

2 

BLM  predicted  =  0.90*ob3erved  -  0.21      r    =  .70, 

2 

modified  P-model      predicted  =  -0.74*observed  -  1.15    r    =  .30. 

The  BLM-predicted  cooling  curve  fit  the  observed  curve  closely,  but 
the  BLM  consistently  predicted  about  3  C  too  warm  for  the  Gainesville 
site.  The  results  of  t-tests  showed  that  the  difference  between  observed 
and  predicted  temperature  was  not  significantly  different  from  zero  for 
the  BLM  at  the  Ruskin  site  and  the  P-model  at  the  Gainesville  site,  but 
significantly  different  from  zero  for  the  BLM  at  the  Gainesville  site  due 
to  a  3.35  C  displacement.  However,  the  variation  in  the  hourly  BLM 
predictions  was  about  half  that  of  the  P-model  during  radiation  nights 
for  both  sites.  The  95$  confidence  intervals  for  the  BLM  at  the  Gainesville 
site  averaged  3.4  C,  for  the  P-model  5.8  C.  The  sensitivity  analysis 
showed  that  the  displacement  was  most  likely  due  to  the  interpolated  input 
files  rather  than  any  particular  input  parameter. 

The  greatest  errors  in  the  overall  BLM  predictions  were  due  to 
advection.  The  advective  terms  in  equation  9  are  not  included  in  the 
model  since  the  model  is  simulating  only  one  point  in  the  horizontal. 
Incorporating  the  BLM  into  a  three-dimensional  regional  or  synoptic  scale 
model  would  improve  its  performance  during  advective  events.  This  could 
possibly  be  done  with  the  soundings  presently  acquired  through  AFOS,  but 
VAS  may  prove  to  be  much  more  valuable  due  to  the  higher  resolution  of 
data. 

The  performance  of  the  Boundary  Layer  Model  shows  improvement  over 
the  P-model.  The  model  in  its  present  form  is  recommended  for  prediction 
of  radiation  night  situations.  Presently,  however,  the  inputs  are  not 
sufficient  for  model  predictions  for  every  desired  location  in  the  state, 
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due  to  the  low  sounding  network  density.  The  Gainesville  sounding 
interpolation  resulted  in  a  consistent  bias.  The  interpolation  procedure 
could  be  used  for  other  sites  after  the  bias  is  determined  and  corrections 
are  included  in  the  prediction.  The  fact  that  all  sounding  sites  in 
Florida  au?e  coastal  is  a  potential  problem.  The  land-sea  interface  produces 
mesoscale  processes  not  associated  with  inland  sites,  such  as  nocturnal 
land  breezes.  VAS  may  prove  valuable  for  the  increase  of  input  data,  both 
in  frequency  and  in  resolution. 

The  Modified  P-model  had  an  advantage  over  the  original  P-model  and 
the  3LM  because  it  could  use  inputs  from  any  location  in  Florida.  However, 
the  Modified  P-model  did  not  make  accurate  site-specific  predictions. 
Use  of  a  single  measurement  as  input  to  the  Modified  P-model  made  the 
model  very  susceptible  to  anomalous  data  and  therefore  less  reliable  than 
the  original  P-model  and  the  BLM.  The  presence  of  clouds  and  moisture 
masked  outgoing  surface  radiation  so  the  temperature  measured  by  the 
satellite  was  not  representative  of  the  actual  surface  temperature. 
Atmospheric  corrections  may  reduce  the  error  caused  by  masking. 

The  Boundary  Layer  Model  was  run  in  real-time  on  some  of  .the  frost 
nights  during  the  Winter  of  1984-85,  but  not  as  an  operational  program. 
The  BLM  could  be  run  in  an  operational  environment  as  an  aid  to  forecasters 
making  frost  night  predictions,  in  the  same  way  that  the  P-model  was 
originally  envisioned.  The  predicted  cooling  curves  could  also  be  sent 
directly  to  users  through  the  Florida  Agricultural  Weather  Network.  The 
predictions  could  be  transferred  as  files,  utilizing  the  graphic  capabilities 
of  user's  microcomputers  to  display  the  predicted  temperature  versus  time 
for  their  area. 
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Models  that  simulate  specific  agricultural  practices,  such  as  orchard 
heating,  sprinkler  irrigation,  and  fog  production,  can  be  incorporated 
into  the  Boundary  Layer  Model  to  study  the  effects  of  these  practices  on 
the  upper  layers  of  the  boundary  layer.  The  BLM  lends  itself  well  to  this 
since  it  has  the  separate  surface  air  layer,  into  which  these  other  models 
could  be  fitted.  Some  examples  of  models  that  could  possibly  be  used  are 
sprinkler  models  (Perry  et  al.,  1982;  Gerber  and  Harrison,  196M),  and 
heater  models  (Welles,  Norman,  and  Martsolf,  1979,  1981). 


APPENDIX  A 
DISCUSSION  OF  INDIVIDUAL  NIGHTS 

Error  in  the  difference  between  BUM-predicted  and  observed  temperatures 
was  consistent  for  the  radiative  frost  nights,  but  not  for  the  advective 
nor  anomalous  nights.  The  following  is  a  discussion  of  processes  that 
were  occurring  and  possible  reasons  why  the  model  did  not  do  as  well  on 
these  particular  nights. 

The  three  advective  nights  at  Gainesville  were  January  5,  21,  and  26 
(Figures  30a,  b,  and  c).  The  predicted  curve  follows  the  observed  curve 
rather  closely  until  the  seventh  predicted  hour  for  the  January  26  case. 
At  this  point  the  winds  shifted  from  southwest  to  northwest  and  increased 
in  speed;  the  increased  turbulence  mixed  the  stratified  atmosphere  causing 
a  temporary  warming  trend.  After  the  tenth  forecast  hour  the  cold  advection 
became  evident  and  the  temperature  dropped  quickly,  well  below  the  predicted 
value. 

The  advection  on  January  21  brought  one  of  the  worst  freezes  In  the 
last  100  years  In  Florida  (Figure  30b).  The  advection  had  already  begun 
when  the  Waycross  sounding  was  made.  The  predicted  cooling  was  similar 
to  the  observed  for  the  early  part  of  the  evening  at  Gainesville.  However, 
the  model  assumed  that  the  outgoing  radiation  flux  would  decrease  as  the 
surface  cooled,  decreasing  the  slope  of  the  cooling  curve  later  In  the 
event.  The  observed  curve  shows  a  steady  drop  In  temperature  through  the 
night  as  the  cold  air  quickly  moved  over  the  area.    The  air  mass  appeared 
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Figure  30.     Cooling  curves  from  three  advective  nights, 
Gainesville  location,  a)  Jan.  5,  1985; 
b)  Jan.  21,  1985;  c)  Jan.  26,  1985. 
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Figure  30 — continued 
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to  be  isothermal  with  height,  unlike    a  radiative  frost  event,  as  shovm 
by  the  fairly  smooth  observed  curve.    Strong  winds  acoorapanled  this  cold 
air  mass,  creating  strong  turbulence,  but  due  to  the  isothermal  conditions 
little  turbulent  exchange  of  heat  occurred. 

The  presence  of  an  upper  air  disturbance  over  the  Gulf  of  Mexico 
creates  the  unusual  temperature  pattern  for  the  night  of  January  13  at 
Gainesville  (Figure  31a).  Clouds  intermittently  passed  over  the  site, 
altering  the  net  radiation.  Also,  the  wind  periodically  increased,  most 
notably  beginning  at  the  fifth  and  eighth  forecast  hours.  The  resulting 
increase  in  mixing  raised  the  1.5  m  air  temperature  during  the  episode  of 
higher  wind.  The  model  has  no  way  of  accounting  for  synoptic  scale 
influences  such  as  this  upper  air  disturbance. 

The  night  of  January  22  was  predicted  (by  forecasters)  to  be  as  cold 
as  the  previous  night,  and  the  BLM  anticipated  that  with  its  prediction 
(Figure  31b).  The  observed  temperature  decreased  rapidly  in  the  early 
evening.  As  the  temperature  began  to  level  off,  the  wind  increased  and 
continued  through  the  entire  evening.  Apparently  a  layer  of  moisture  also 
moved  into  the  area  cutting  down  on  the  radiative  heat  loss. 

The  air  temperature  decreased  quite  rapidly  early  on  the  night  of 
January  23  (Figure  31c),  much  more  quickly  than  the  model  predicted. 
However,  the  cooling  was  sharply  cut  off  with  an  increase  in  the  wind  and 
mixing  at  the  fifth  prediction  hour.  The  predicted  minimum  and  the  observed 
temperature  were  within  0.5  C  and  30  minutes  of  each  other  at  the  end  of 
the  frost  event. 

For  seven  of  the  eleven  nights  studied  at  the  Huskin  site,  the  observed 
curve  followed  the  typical  decaying  exponential  temperature  drop  of  a 
radiative  frost  night.    However,  unlike  the  Gainesville  site  which  tends 
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Figure  31.     Cooling  curves  from  anomalous  nights,  Gaines- 
ville location,     a)  Jan.  13,   1985;  b)  Jan.  22, 
1985;  c)  Jan.  23,  1985. 


93 


c 


Figure  31 — continued 
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to  be  a  cold  pocket  on  frost  nights,  the  Huskin  site  is  located  near  the 
Tampa  Bay.  The  influence  of  this  large  body  of  water  can  be  seen  in  the 
observed  curves  of  the  graphs.  The  prevailing  winds  carried  some  moisture 
off  the  bay  over  the  site  in  some  cases,  partially  masking  the  outgoing 
radiation.  The  site  is  in  close  enough  proximity  to  the  bay  to  be  influenced 
by  the  nocturnal  land  breeze  that  develops  because  of  the  land-sea 
temperature  differential.  This  breeze  could  have  been  responsible  for  an 
increase  in  mixing  of  the  stratified  layers.  The  radiation  and  the  land 
breeze  were  probably  responsible  for  the  environment  not  cooling  quite  to 
the  degree  that  the  BLM  predicted.  This  occurred  on  the  nights  of  December 
8,  1934,  January  23,  1985,  and  February  15,  1985  (Figures  32a,  b,  and  c, 
respectively).  The  night  of  January  23  is  a  good  example  of  possible  bay 
influence.  The  BLM  predicted  strong  radiative  heat  loss  in  the  early 
hours,  then  a  rapid  decrease  in  heat  loss  after  the  second  forecast  hour. 
The  actual  cooling  was  moderated  in  the  early  hours,  but  then  Increased 
after  the  seventh  forecast  hour  (Figure  32b).  The  predicted  minimum  was 
only  about  1  C  less  than  the  actual  minimum. 

The  BLM  predicted  the  cooling  quite  accurately  for  the  nights  of 
December  9,  1984,  January  24,  1985,  and  January  27,  1985  (Figures  33a,  b, 
and  c,  respectively).  It  is  interesting  to  note  how  close  the  predicted 
and  observed  curves  fit,  particularly  on  January  24.  In  this  case  the 
model  predicted  the  strong  radiative  heat  loss  early  in  the  evening  that 
is  quite  evident  on  the  observed  curve.  The  inflection  point,  where  the 
strong  cooling  decreases  considerably,  differed  by  less  than  one  hour 
between  the  two  curves. 

Three  of  the  ten  nights  at  the  Ruskin  were  affected  by  cold  advection, 
December  7,  1984,  January  21,  1985,  and  January  26,  I985  (Figures  34a,  b. 
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Figure  32.     Cooling  curves  from  anomalous  nights,  Ruskin 
location,     a)  Dec.  8,  1984;  b)  Jan.  23,  1985; 
c)  Feb.  16,  1985. 
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Figure  32 — continued 
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Figure  33.     Cooling  curves  from  radiation  nights,  Ruskin 
location,     a)  Dec.  9,  1984;  b)  Jan.  24,  1985; 
c)  Jan.  27,  1985. 
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Figure  33 — continued 
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Figure  34.  Cooling  curves  from  advection  nights,  Ruskln 
location,  a)  Dec.  7,  1984;  b)  Jan.  21,  1985 
c)  Jan.  26,  1985. 


100 


RUSK IN  626 


sxssstxxtsssxtx 

TINC  (NOUKS) 

c 


Figure  34 — contlaued 
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and  c).  On  December  7,  the  predicted  curve  was  similar  to  the  observed 
curve  until  the  fifth  prediction  hour,  then  the  cold  air  mass  began  moving 
in  and  the  temperatures  dropped  rapidly.  On  January  21  the  colder  air 
arrived  just  after  the  sounding  was  made.  This  proved  to  be  the  worst 
3LM  forecast,  because  the  unaffected  air  mass  over  Ruskin  at  the  sounding 
time  was  rather  humid,  and  the  surface  was  warm.  The  BLM  anticipated  a 
very  slow  drop  in  temperatures  through  the  night,  when  in  actuality  the 
advected  air  dropped  the  temperatures  quickly.  The  prediction  for  Jamuary 
26  was  quite  close  to  the  observed  up  until  the  eleventh  forecast  hour 
when  the  advection  began  to  occur. 


APPENDIX  B 
ADDITIONAL  GRAPHS  OF  COOLING  CURVES 
FOR  INDIVIDUAL  NIGHTS 
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APPENDIX  C 

LISTING  OF  PROGRAMS  AND  PROCEDURE  TO  PRODUCE 
TEMPERATURE  PREDICTIONS  FROM  THE  MODELS 

1 .    The  P-model 

Equipment  used:    HP  21MX-E 

HP  terminal 
Programs  used:  AKD 

PMODG  (version  of  P-model  used  in  this  study) 

The  P-model  input  files  must  be  in  a  special  non-ASCII  form  to  be 
read  by  the  model  software.  The  measurements  can  be  placed  in  an  ASCII 
file  in  the  following  format: 


year, Julian  day 
station 

time,T3ur,Ts10,T350,Ta1 .5,Ta3,Ta9,v,dir,Rn 

where  Tsur  =  soil  surface  temperature  (F), 
TslO  =  10  cm  soil  depth  temperature, 
T35O  =  50  cm  soil  depth  temperature, 
Ta1.5  =  1.5  m  air  temperature, 
Ta3  =  3«0  m  air  temperature, 
Ta9  =  9.0  m  air  temperature, 
V  =  wind  speed  (mph), 
dir  =  wind  direction 
Rn  =  net  radiation  (Cal  cm~    s"  ). 

example: 

013,1985 
GNV 

1,46.4,57.6,53.5,40.3,41.0,41.5,7.6,0,1.0 
2,45.1,57.6,52.7,38.8,39.4,39.9,5.7,0,1.0 
3,43.3,57.6,52.0,35.6,37.0,37.9,5.7,0,1.0 

Running  program  AKD  will  convert  this  ASCII  file  to  the  P-model  input 

file.    The  file  name  usually  takes  the  form  KYYDDD,  where  yy  is  the  year 

and  ddd  the  Julian  day  (e.g.  K85OI3).    When  the  P-model  is  run  interactively. 
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the  user  will  be  interrogated  for  the  input  file  name.  At  the  present 
time  AKD  must  have  data  from  all  10  stations.  However,  the  desired  station 
is  the  only  one  that  needs  real  data,  and  the  user  will  be  interrogated 
for  the  number  of  that  station. 

2.    The  Modified  P-model 

Equipment  used:      HP  21  MX-E 

HP  terminal 
Programs  used:       GOES,  SCTRZ 

EDGE 

LOCAT  (or  LOG) 

SAT2B 

AKD 

PMODG 

The  programs  GOES  and  SCTRZ  are  scheduled  by  the  SFFS  scheduler  to 
produce  the  Florida  sector.  EDGE  and  LOCAT  (or  LOG)  are  programs  used 
for  registration  of  the  satellite  maps.  EDGE  enhances  the  strongest 
temperature  gradients  so  the  outline  of  Lake  Okeechobee  stands  out.  LOCAT 
compares  the  Lake  outline  with  a  base  map  to  determine  the  offsets.  LOG 
matches  the  Florida  sector  grid  outline  with  a  base  map,  to  be  used  if 
Lake  Okeechobee  is  covered  with  clouds.  All  three  programs  are  run 
interactively  and  will  ask  for  the  map  names. 

SAT2B  reads  the  satellite  temperatures  and  produces  an  ASCII  file 
for  AKD  to  process  into  P-model  input  files.  The  user  must  create  a  file 
with  the  map  names  and  x  and  y  offsets  for  SAT2B  to  read.  SAT2B  will  ask 
for  the  file  name.    The  procedure  from  this  point  on  is  the  same  as  in  1. 


3.    The  Boundary  Layer  Model 

Equipment  used:    HP  21  MX-E 
HP  terminal 
CYBER  730  or  VAX  750 
VIKING,  VT,  or  HP  terminal 

Programs  used:  DSORT 
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UPAIR 

IN2 

EMUL8 

PBL  (latest  version  of  BLM  on  CYBER)  or 
BLACKS  (latest  version  of  BLM  on  VAX) 

Upper  air  soundings  are  acquired  from  Ruskin  WSO  through  the  DS-1000 

link  by  program  DSORT.    They  may  need  minor  editing,  and  their  format 

should  be  identical  to  the  exsunple  in  Table  3.    All  seven  soundings  should 

be  obtained.    Program  UPAIR  reads  the  upper  air  soundings,  strips  off  the 

necessary  values,  sorts  the  values,  and  places  them  into  files.  Program 

IN2  then  reads  those  files,  interpolates  the  data,  and  produces  BLM  input 

files  for  each  station,  including  the  additional  inputs  found  in  Table  2. 

These  programs  are  run  on  the  HP  21  MX-E.    The  files  need  to  be  transferred 

to  the  CYBER  or  the  VAX  where  the  executable  BLM  program  is  located.  This 

can  easily  be  done  with  the  HP  program  EMUL8,  through  file  transfer.  Once 

the  input  files  are  transferred,  the  best  approach  is  to  use  a  VIKING 

terminal  for  the  CYBER  or  a  VT  terminal  for  the  VAX  to  run  the  BLM.  The 

BLM  is  interactive,  and  will  interrogate  the  user  for  the  desired  input 

file. 


APPENDIX  D 
BLM  AND  P-MODEL  INPUT  FILES 
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P-model  input  files 


ONV 

005, 1985 

1  ,51  .  1  ,60. 6,65. 6, 115.7,  *6,0,*6. 2,6. 0,0,  1 .0 
2, 50. 5, 59.1, 65. 1,**.  4, 45. 0,115.0, 1.0, 0,1.0 
3. 19. 1,58. 5, 6*. 9, 43.3.43. 9, 41. 1,3. 2, 0.1.0 

OHT 

006, 1985 

1 ,49.3,56.1,60.8,42.4,44.4,45.7,3.0,0, 1.0 
2,46.4,54.9,60.6,  34.7,  35.8,36.7,  1  .0,0,  1  .0 
3, 44. 2, 53. 6, 60. 6, 31. 5, 34. 0,34. 5, 0.0, 0,1.0 

ONT 

013,1985 

1  , 46. 4, 53. 6, 57. 6, 40. 8, 41  . 0,41  . 5, 4. 0,0,  1  .0 
2,45.  1  ,52.7,57.6,  38.8,39.4,39.9,  3.0,0,  1  .0 
3, 43. 3, 52. 8, 57. 6,  35. 6,  37. 0,37. 9, 3. 0,0,  1  .0 

ONT 

016, 1985 

1  ,50.0,56.8,55.9.45.  1  ,47.5,49.  1  ,0.0,0,  1  .0 
2,45.7,55.4,56. 1 ,37.6,40.6.40.8.0.0,0, 1 .0 
3,43.5.54.1,55.9,33.8.37.4,38. 1,0.0,0, 1.0 

GMT 

021 ,1985 

1  ,43.3,50.9,55.4,36.5,36.5,  36.5,  18.0,0,  1  .0 
2.40.8,49.6,55.2,32.0,32.0,32.0,24. 1 ,0, 1 .0 
3,38.7.48.4,55.2,29.5,29.5,29.5,22.5,0, 1 .0 

GMT 

022, 1985 

1  ,35.  1  ,42.8.51  .  1  .25.2,25.9,26.1  ,8.3.0,  1  .0 
2, 33. 4. 41. 9. 51. 3,  18. 9, 21. 9, 23. 5. 0.0, 0.  1  .0 
3, 32. 6,  40.  1,41. 1  ,14. 2, 15. 6, 19. 2, 0.0. 0,1.0 

OUT 

023,1985 

1  , 42. 8, 45. 9, 48. 6. 41  . 9. 42. 6, 43. 5, 7. 8. 0,  1  .0 
2.39.7.45.  1  ,48.7,  38.5.39.4,39.7.8.8,0,  1  .0 
3,37.4,44.1,48.9,35.6,36.5,36.9.12.0,0,1.0 

ONT 

024, 1985 

1  ,47.3,50.0,49.  1  ,47.8,48.7,  49.5.9.0,0,  1  .0 
2,43.7,49.1,49.3,44.1,45.3,45.9,8.5,0,1.0 
3, 41. 0,47. 7,  49. 3, 41. 0,42. 4, 43.  3, 4. 0,0, 1.0 


OHT 

026, 1985 

1,54.5,56.8,52.5,56.5,57.0,57.6,  10.0,0,  1  .0 
2,51  .6,55.6,52.7,53.8,54.5,54.9,9.0,0,  1  .0 
3, 50. 5, 54. 3, 52. 9, 52. 7, 53. 4, 54. 0,12. 0,0, 1.0 

ONV 

027, 1985 

1,48.2,53.4,52.5.41.2,43.0,44.2,0.0,0,1.0 
2,43.0,51.8,52.7,31.3,34.5,35.4,0.0,0,1.0 
3,39.9,50.2,52.9,28.2,30.7,30.7,0.0,0, 1.0 

GNT 

041 , 1985 

1,51.8,56.8,56.7,49.3,49.6,50.7,11.4,0,1.0 
2,46.8.55.4,56.7,37.2,40.6,41.4,0.0,0,  1  .0 
3,43.5,53.6,56.8,32.7,33.3,34.9,0.0,0,1.0 

GNT 

048, 1985 

1.54.1,57.5.52.9.51.1,52.0,53.2,0.0,0.1.0 
2. 48. 6, 56. 5. 53. 1.42. 1,45. 3, 46. 9, 0.0, 0,1.0 
3.45.3.54.9.53.1,39.2,41.2,43.0,0,0,0,1.0 
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